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DESIGN  CONSIDERATIONS  FOR  AN  EXPERIMENTAL 
HEATED  FLAT  PLATE  STABILITY  STUDY 


Abstract 
by 
ROBERT  LEROY  LOWELL,  JR. 

The  design  and  construction  of  an  apparatus  for  experimentally 
examining  the  stability  of  a  uniform,  parallel  water  flow  over  a 
heated,  constant-temperature,  sharp-edged,  flat  plate  to  a 
mechanically  generated  disturbance  is  discussed.  Consideration 
is  given  to  the  effects  of  the  streamwise  pressure  gradient,  test 
section  turbulence  intensity  level,  uniformity  of  the  flow  and 
orientation  of  the  plate  in  it,  material  considerations  for  the 
plate,  maximum  utilization  of  the  existing  facility,  heat  genera- 
tion required  to  maintain  a  constant  temperature  plate  for  different 
flow  speeds,  and  production  of  a  small  disturbance  of  specified 
amplitude.  The  last  two  aspects  of  the  analysis  are  mathematically 
modelled  and  the  resulting  analytical  solutions  used  to  establish 
the  strength  and  distribution  of  heat  sources,  for  the  former,  and 
the  effect  of  viscous  damping,  resonance  frequencies,  tension,  etc., 
on  a  vibrating  ribbon,  for  the  latter.  Where  feasible,  experimental 
checks  are  made  to  verify  that  the  designed  plate  assembly  operates 
as  required. 
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CHAPTER  1 
INTRODUCTION 

Since  the  beginning  of  the  twentieth  century,  much  interest 
has  been  expressed  and  effort  expended  in  the  study  of  conditions 
affecting  the  stability  of  laminar  boundary  layers  as  a  possible 
explanation  for  the  onset  of  boundary  layer  turbulence  (see  refer- 
ence (1)  for  a  comprehensive  bibliography).   One  such  condition  - 
the  possible  stabilizing  effect  of  heating  due  to  a  liquid's  vis- 
cosity-temperature dependence  -  provides  the  motivation  for  this 
investigation.   Choosing  a  relatively  simple,  well-defined  experi- 
mental configuration  -  that  of  a  uniform,  parallel,  water  flow 
over  a  heated,  constant-temperature,  sharp-edged,  flat  plate  - 
the  significance  of  this  condition  is  to  be  examined.   It  is  anti- 
cipated that  the  empirical  results  obtained  from  experiments  using 
the  apparatus  described  below  might  be  reasonably  compared  with  the 
analytical  solution  obtained  by  Wazzan,  Okamura,  and  Smith  (2,3). 

The  ensuing  discussion,  then,  outlines  the  various  facets  of 
design,  construction,  and  where  feasible,  testing  involved  in  the 
model  simulation.   Since  a  stability  study  is  to  be  conducted,  not 
only  should  the  conditions  delineated  above  for  the  model  be  dupli- 
cated as  closely  as  possible,  but  all  those  masking  the  effect  of 
the  heated  flat  plate  must  be  reduced  to  a  level  of  insignificance. 


For  the  latter,  such  conditions  included  the  effects  of  surface 
roughness  and  waviness  of  the  plate,  the  effect  of  a  streamwise 
pressure  gradient,  and  the  possible  importance  of  buoyancy  effects. 
Eventually  each  of  these  influenced  design  considerations,  as  did 
the  more  practical  requirements  concerning  the  selection  of  the 
proper  materials,  ease  of  access  to  the  parts  of  the  final  assembly, 
replaceability  of  components,  elimination  of  any  extraneous  vibra- 
tions, maximum  utilization  of  the  available  testing  length,  and  the 
paramount  requirement  that  the  assembly  be  able  to  seal  against 
water  leakage. 

Specifically,  the  design  aspects  for  the  experimental  model 
were  allocated  to  three  basic  areas: 

1.  the  simulation  of  a  flat  plate  at  zero  incidence  to  a 
uniform,  parallel  flow  with  very  low  turbulence  intensity  and  zero 
pressure  gradient; 

2.  the  construction  of  a  flat  plate  that  would  give  constant 
temperature  on  the  working  surface  (the  portion  of  the  plate  adja- 
cent to  the  flow  over  which  the  boundary  layer  will  build  and 
measurements  will  be  taken)  with  convective  heat  transfer  along  its 
length; 

3.  the  design  of  a  device  for  introducing  a  disturbance  into 
the  boundary  layer  of  the  plate.   Note  that  no  restriction  is  placed 
on  the  amplitude  of  the  disturbance,  so  that  both  finite  oscillat- 
ing waves  and  the  "infinitesimal  Tollmein-Schlichting"  ones  can  be 
generated. 


By  presenting  several  possible  methods  for  simulating  a 
given  condition,  and  then  discussing  the  criteria  governing 
the  selection  of  a  particular  one,  insight  is  given  to  the 
interdependency  of  the  numerous  design  aspects.   In  cases 
where  discrepancies  inherently  exist  between  a  proposed  com- 
ponent design  and  the  assumed  experimental  configuration, 
modifications  to  the  former  are  made  to  optimize  overall 
system  performance.   Since  the  systematic  search  for  per- 
formance requirements  for  the  composite  apparatus  followed  a 
logical  progression  once  all  the  required  features  were 
identified,  the  ensuing  discussion  will  follow  that  same  order, 


CHAPTER  2 

SIMULATION  OF  A  FLAT  PLATE  AT  ZERO  INCIDENCE  TO  A 
UNIFORM,  PARALLEL  FLOW  WITH  LOW  TURBULENCE  LEVEL 
AND  ZERO  PRESSURE  GRADIENT 


2.1  Testing  Facility 

A.   Physical  characteristics  and  operating  capabilities 

The  experimental  stability  investigation  outlined  in  Chapter  1 
is  to  be  conducted  in  the  6"  x  9",  rigid-walled,  constant-area 
test  section  of  the  closed-loop,  low-turbulence  water  tunnel  at 
Case  Western  Reserve  University.   Designed  and  constructed  by 
Domholdt  (4),  with  a  25-to-l  contraction  ratio  nozzle  located 
downstream  of  a  series  of  damping  screens,  the  facility  is  capable 
of  operating  at  test  section  free  stream  speeds  of  up  to  forty 
feet  per  second  with  longitudinal  turbulence  intensities  levels 
ranging  from  .04%  at  5  fps  to  .15%  at  35  fps.   With  a  utilizable 

testing  length  of  12t'\  this  reflects  a  maximum   Ke<-*  =  3147. 

3 

The  anodized  aluminum  flat  plate  (9"  wide,  ■$     thick,  and 

o 

12t-"  long)  over  which  measurements  are  to  be  taken  is  fitted  into 

a  plexiglass  frame  mounted  on  top  of  the  test  section  so  that  its 

9 
sharp  leading  edge  extends  into  the  flow  approximately  ~ . 

lo 

Figure  2.1  shows  the  location  of  the  plate  and  the  portion  of  the 
tunnel  circuit  germane  to  this  discussion.   Further  explanation 
for  the  choice  of  this  particular  experimental  configuration  con- 
stitutes the  remainder  of  this  paper. 
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B.   Uniformity  and  turbulence  intensity  of  the  flow 

Before  proceeding  with  a  discussion  of  the  plate  itself,  it 
would  seem  feasible  to  examine  briefly  the  attainability  and 
ramifications  of  the  requirement  for  a  uniform,  parallel  flow  with 
a  low  turbulence  intensity. 

Due  to  the  symmetry  of  the  nozzle  section  immediately  preceding 
the  test  section  in  the  water  tunnel,  and  the  fact  that  the  leading 
edge  of  the  flat  plate  would  be  located  approximately  1— _"  downstream 
of  the  nozzle-test  section  juncture,  it  is  expected  that  any  second- 
ary "cross"  flows  arising  from  the  curvature  of  the  nozzle  would  be 
insignificant,  and  the  flow  would  approach  the  location  of  the 
leading  edge  in  a  uniform,  parallel  fashion.   This  expectation  is 
confirmed  by  Domholdt  (4)  and  Nice  (5),  whose  pitot  tube  and  ane- 
mometry  equipment  measurements,  respectively,  verified  an  essentially 
slug-like  flow  across  the  test  section  for  all  speeds,  with  the 
expected  boundary  layer  retardation  at  the  sides.   The  former  reports 
a  uniformity  of  to  within  +  0.5%  over  the  entire  speed  range,  but 
notes  that  in  most  cases,  better  than  +  0.25%  was  observed  (see 
figures  6-10  of  reference  (A)). 

Concerning  the  free  stream  turbulence  intensity,  Schlichting  (6) 
points  out  that  based  on  the  results  of  Schubauer  and  Skramstad  (7), 
unless  the  turbulence  level  of  the  flow  is  sufficiently  small,  transi- 
tion will  not  be  caused  by  the  amplification  of  artificial  oscilla- 
tions of  certain  frequencies,  but  by  the  random  disturbances  already 


present  in  the  flow.   Specifically,  for  a  turbulence  intensity 
above  .03%,  transition  is  caused  directly  by  random  disturbances 
rather  than  by  the  selective  amplification  of  sinusoidal  oscilla- 
tions.  If  on  the  other  hand,  the  flow  turbulence  intensity  is  low, 
then  amplified  artificially  introduced  sinusoidal  disturbances  of 
a  definite  frequency  can  be  detected  even  in  the  presence  of  natural 
oscillations  inherent  in  the  flow.   Schlichting  (6),  using  the 
results  of  Granville  (8) ,  has  shown  that  by  combining  the  experi- 
mental results  of  Schubauer  and  Skramstad,  Hall  and  Hislop,  and 
Dryden,  the  difference  between  the  Reynolds  numbers  at  the  points 
of  transition  and  instability  (plotted  as  a  function  of  the  tur- 
bulence intensity)  decreases  with  increased  turbulence,  but  does 
not  equal  zero  until  the  turbulence  intensity  is  approximately  2% 
to  3%. 

Since  the  establishment  of  low  turbulence  levels  is  so 
important  for  stability  measurements,  an  attempt  to  reduce  the  level 
of  the  testing  facility  still  further  has  been  made  by  the  addi- 
tion of  four  more  damping  screens  to  the  original  seven,  but 
little  or  no  decrease  was  observed  (5).   Apparently,  as  discussed 
by  Schubauer  and  Skramstad  (7),  the  noise  generated  as  a  result 
of  the  increased  power  needed  to  force  the  flow  through  the  additional 
screens  at  the  same  flow  rate  offset  any  reduction  in  free  stream 
turbulence  intensity.   It  should  be  noted  that  the  turbulence  inten- 
sity measurements  of  Domholdt  (A)  and  Nice  (5)  are  only  one- 
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dimensional  ones  made  in  the  direction  of  greatest  fluctuation 
(longitudinally)  and  then  extrapolated  to  three-dimensional,  and  as 
such,  should  be  considered  as  upper  bounds  on  the  actual  free  stream 
turbulence  intensities. 

2.2  Orientation  of  the  Plate  in  the  Test  Section 
A.   General  location  of  the  plate 

Now  that  it  is  determined  that  a  uniform,  parallel,  low- 
turbulence  intensity  flow  can  be  established  in  the  test  section 
(section  2. IB),  the  simulation  of  a  flat  plate  at  zero  incidence 
to  this  flow  is  examined.   Logically,  the  first  problem  encountered 
concerned  the  general  orientation  of  the  flat  plate  in  the  test  sec- 
tion, for  this  ultimately  must  determine  the  plate  design. 

In  an  attempt  to  utilize  the  existing  test  section  -  one  in 
which  no  provision  had  been  built  to  compensate  for  the  favorable 
pressure  gradient  that  would  result  from  boundary  layer  build-up 
(e.g.,  moveable  test  section  walls  to  change  the  downstream  cross- 
sectional  area,  boundary  layer  suction  along  the  walls,  etc.)  - 
the  plate  was  not  put  in  the  center  of  the  test  section,  contrary 
to  the  majority  of  wind  tunnel  experiments.   To  do  so  would  cause 
the  boundary  layer  growth  on  both  sides  of  the  plate  as  well  as 
the  test  section  walls,  thereby  decreasing  the  already  limited  effective 
potential  flow  area  and  so  increasing  even  more  the  streamwise 
pressure  gradient.   As  already  noted,  most  of  existing  experimenta- 
tion for  the  type  problem  considered  here  has  been  carried  out  in 


a  gaseous  medium  using  a  wind  tunnel.   For  a  water  tunnel  where 
sealing  and  leakage  problems  arise,  as  much  as  possible  of  the 
associated  equipment  for  the  plate  normally  placed  in  the  tunnel 
on  the  non-testing  side  of  the  plate  now  must  be  on  the  outside. 
Therefore,  it  is  felt  that  ease  of  direct  access  to  at  least 
one  side  of  the  plate  would  be  advantageous. 

A  second  possibility  for  the  placement  was  to  put  the 
plate  on  the  bottom  of  the  test  section.   However,  since  it  was 
anticipated  that  the  water  tunnel  should  at  times  be  operated 
at  low  flow  rates,  and  that  stability  measurements  would  be 
taken  over  the  heated  flat  plate,  the  Rayleigh  instability 
problems  (where  a  top-heavy  fluid  is  statically  stable  under  the 
joint  influence  of  conduction  and  viscosity  until  the  vertical 
temperature  gradient  becomes  sufficiently  large)  normally  associated 
with  free  convection  flow  might  become  significant.   Of  course, 
similar  reasoning  followed  with  respect  to  locating  the  flat  plate 
on  either  of  the  side  walls. 

Therefore,  the  flat  plate  was  installed  on  the  top  of  the 
test  section.   In  this  way,  external  access  was  granted  to  the 
heating  elements,  pressure  taps,  thermocouples  imbedded  in  the 
plate,  and  parts  of  the  perturbing  mechanism. 
B .   Amount  of  plate  insertion  into  the  test  section 

This  choice  of  position  for  installation,  however,  created 
other  problems  when  trying  to  simulate  the  actual  flow  phenomenon 
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used  for  the  model  -  problems  not  present  in  most  wind  tunnel 
experiments  as  already  noted.   For  example,  in  the  latter  facility 
where  a  plate  could  be  placed  so  that  free-stream  parallel  streamlines 
would  split  the  leading  edge,  and  test  section  walls  could  be 
adjusted  to  accommodate  the  pressure  gradient,  the  flow  could  be 
channelled  around  the  plate  so  as  to  simulate  a  flat  plate  at 
zero  incidence  to  a  uniform,  parallel  flow. 

On  the  other  hand,  for  this  analysis,  so  that  the  boundary 
layer  developing  in  the  nozzle  section  of  the  water  tunnel  did 
not  interfere  with  that  on  the  plate,  the  requirement  was  im- 
posed that  the  plate  be  inserted  into  the  test  section  to  a 
distance  greater  than  the  anticipated  nozzle-section  boundary 
layer  thickness.   The  undetermined  quantity,  then,  was  how  great 
this  distance  should  be  to  ensure  that  the  leading  edge  was  at 
least  into  the  region  of  uniform  or  "slug-like"  flow.   According 
to  Domholdt  (4),  although  the  boundary  layer  thickness  (based  on 
a  2%  velocity  defect)  varied  considerably  at  test  section  speeds 
less  than  15  feet  per  second  (fps),  it  remained  essentially  con- 
stant (and  greater  than  for  the  lower  speed  range)  at  about  0.14 
inches  for  test  section  speeds  above  15  fps  (see  figure  11, 
reference  (4)).   His  explanation  for  this  observation  was  that 
the  boundary  layer  remained  laminar  up  to  a  speed  of  at  least 
7  fps  (where  £  °C  U^  );  but  beyond  about  15  fps,  the  point  of 
transition  apparently  moved  farther  upstream  than  it  would  have 
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for  a  flat  plate  for  the  same  velocity  increase.   Consequently,  the 
boundary  layer  thickness  remained  essentially  constant  instead  of 
thinning  as  it  would  for  a  fixed  position  on  a  flat  plate  (fioC  U^  ) 
Therefore,  in  a  compromise  between  assuring  that  all  the  pre- 
existing boundary  layer  could  be  removed  at  the  leading  edge  of 
the  plate  (in  a  manner  discussed  in  the  next  section),  and  the  de- 
sire to  remove  the  least  amount  of  water  at  the  leading  edge,  the 
plate  was  initially  extended  into  the  test  section  3/16"  or 
approximately  .048"  greater  than  the  maximum  expected  nozzle-exit, 
boundary  layer  thickness.   It  should  be  noted  that  the  design  of 
the  plate  assembly  allowed  for  increasing  this  distance  by  inser- 
tion of  different  size  gaskets  between  the  plate  and  plexiglass 
frame  holding  it  in  place  (see  figure  2.2),  if  later  testing  deemed 
it  necessary.   Eventually  this  distance  was  increased  to  approxi- 
mately 9/16"  in  an  attempt  to  reduce  the  test  section  streamwise 
pressure  gradient,  as  discussed  in  section  2.4C. 

2.3  Leading  Edge  Suction  Requirements 
A.   Pumping  and  piping  considerations 

With  the  irregular  constriction  now  presented  to  the  test  sec- 
tion flow  by  the  protruding  flat  plate,  some  provision  had  to  be 
made  to  ensure  that  the  flow  remained  uniform  and  parallel  as 
before.   The  obvious  solution  was  to  provide  just  the  right  amount 
of  suction  at  the  leading  edge  of  the  plate  so  that  all  the  flow 
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between  the  bottom  of  the  plate  and  the  top  of  the  test  section 
(Including  especially  the  boundary  layer  build-up  from  the  nozzle) 
would  be  channelled  off,  thereby  permitting  only  uniform,  parallel 
flow  on  the  testing  side  of  the  plate  at  the  leading  edge.   Hence, 
streamlines  in  the  plane  of  the  plate  would  split  its  leading  edge 
and  a  Blasius  type  boundary  layer  would  form.   Note  that  "suction" 
is  used  here  to  denote  only  the  process  for  ensuring  that  all  the 
flow  approaching  the  plate  between  the  top  of  the  test  section  and 
the  leading  edge  be  channelled  into  the  "suction  channel,"  and 
should  not  be  confused  with  the  term  "boundary  layer  suction"  in  its 
normal  context.   By  applying  suction  with  the  protruding  plate,  it 
is  anticipated  that  any  irregularities  in  the  boundary  layer 
(e.g.,  caused  by  the  discontinuous  nozzle-test  section  -juncture 
or  a  dye  injection  device  to  be  discussed  in  the  next  section) 
can  be  drawn  off  also. 

Similar  attempts  at  this  type  of  flat  plate  design  (9,10) 
have  resulted  in  a  limitation  on  the  maximum  attainable  leading 
edge  suction,  thereby  restricting  useable  test  section  operating 
speed  to  less  than  ten  f ps ,  or  about  one-quarter  of  the  maximum 
capacity  of  the  tunnel.  It  was  felt  that  this  limitation  was  a 
direct  result  of  inadequate  pumping  facilities  for  handling  the 
volume  of  flow  required  at  the  net  positive  suction  head  (NPSH) 
that  existed  at  the  greater  speeds.   Perhaps  this  particular 
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aspect  of  the  design  required  more  attention  than  might  otherwise 
have  been  assumed. 

The  whole  design  and  construction  of  the  leading-edge  suction 
apparatus  was  based  upon  making  the  head  loss  through  that  system 
as  small  as  possible  in  order  to  be  able  to  use  the  smallest  (and 
consequently,  the  least  expensive)  pump.   Since  the  large  volume 
of  water  that  had  to  be  withdrawn  made  it  impractical  to  discard 
the  deionized  water  of  the  main  tunnel  after  drawing  it  off,  it 
was  decided  that  this  water  should  be  reinserted  into  the  main 
tunnel  loop.   In  this  way,  the  suction  apparatus  would  just  form 
an  auxiliary  loop  of  the  main  tunnel  flow.   It  was  required  that 
the  water  be  reinjected  in  the  tunnel  at  a  position  where  the 
static  pressure  was  relatively  low  and  where  the  injected  flow  would 
not  affect  the  flow  through  the  test  section.   Examination  of 
figure  13  of  reference  (4)  revealed  that  the  lowest  static  pressure 
in  the  tunnel  during  operation,  excluding  the  test  section  itself, 
was  in  the  transition  section  just  upstream  of  the  main  pump. 

It  was  anticipated  that  this  auxiliary  loop  should  have 
between  the  test  and  transition  sections  a  pump  to  circulate  the 
water  (essentially  just  overcoming  the  head  losses  in  the  system), 
one  or  more  valves  to  regulate  the  flow  through  the  piping,  and 
some  sort  of  device  to  measure  the  volume  of  water  which  would 
cover  the  full  range  of  flows  needed.   Here  again,  the  selection 
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of  the  valves  and  piping  size  represented  a  compromise  between 
the  higher  head  losses  at  large  flow  rates  through  smaller  piping, 
and  higher  costs  and  installation  inconvenience  of  larger  size 
piping  and  fittings.   After  some  calculation,  it  was  determined 
that  a  pipe  size  of  2-1/2"  with  associated  fittings  and  the 
physical  location  of  each  component  would  result  in  an  NPSH  of 
about  8.7'  and  require  a  pump  of  at  least  4  BHP.   With  these 
figures  in  mind,  it  became  obvious  why  difficulty  had  previously 
been  experienced  in  obtaining  sufficient  suction  for  tunnel 
operation  at  speeds  above  ten  fps  (where  two  3/4  HP  pumps  were 
used  in  parallel).   The  size  of  the  pump  was  chosen  large  enough 
to  allow  for  the  possible  extension  of  the  plate  further  into 
the  test  section  (and  hence  greater  flow  rate  pumping  requirement) 
for  reasons  discussed  previously.   Consequently,  a  Goulds 
Model  3643,  2-1/2  x  2-1/2,  close  coupled,  single  stage,  enclosed 
impeller,  centrifugal  pump  rated  at  5  HP  at  3500  RPM  was  chosen. 
It  should  be  noted  that  a  centrifugal  pump  was  chosen  rather  than 
other  types  for  the  following  reasons:   quiet  operation;  relatively 
low  initial  cost  and  inexpensive  and  easy  maintenance;  small  size 
for  associated  output;  no  excessive  pressures  even  with  the 
discharge  valve  closed;  smooth  non-pulsating  flow;  simple  construc- 
tion with  relatively  few  moving  parts  and  therefore  less  chance 
of  failure;  and  high  efficiency  for  relatively  large  volumes  of 
liquid  at  small  heads. 
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Since  the  size  of  the  valve  required  for  this  size  piping  would 
hardly  provide  the  fine  adjustment  needed  to  regulate  the  amount  of 
suction  needed  at  the  leading  edge  of  the  plate,  it  was  decided  that 
a  fine  control  valve  would  be  needed,  so  a  1-1/2"  needle  valve  was 
placed  in  parallel  with  a  3"  globe  valve. 

The  final  design  of  the  system  appears  schematically  in 
figure  2.1.   All  piping  with  the  exception  of  that  on  the  branch  with 
the  needle  valve  is  of  2-1/2"  polyvinyl  chloride  (PVC)  (the  latter 
being  1-1/2"  PVC  piping).   Additional  considerations  in  the  place- 
ment of  valves  included  the  requirements  that  the  system  not  have 
to  be  completely  dismantled  when  the  water  tunnel  has  to  be  partially 
disassembled  (to  check  the  screens,  remove  the  test  section,  etc.), 
couplings  be  strategically  placed  so  as  to  allow  easy  maintenance 
and  disassembly  when  necessary,  the  minimum  amount  of  piping  be  used 
to  reduce  the  head  loss,  and  the  valves  and  rotameter  be  placed 
so  as  to  be  easily  accessible  from  one  central  location  -  the 
operator's  stand. 
B.   Attainment  of  the  proper  amount  of  suction 

Due  to  the  wide  variation  of  flows  that  would  be  measured  for 
the  full  range  of  tunnel  operation,  it  was  necessary  to  devise  a 
means  for  measuring  approximately  how  much  liquid  was  being  pumped 
through  the  suction  system.   Because  the  flow  in  the  tunnel  could 
be  measured  but  not  exactly  reproduced  from  one  run  to  the  next, 
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the  exact  suction  needed  for  any  given  experiment  could  not  be 
obtained  by  lust  reading  a  dial  or  gauge  in  the  suction  line. 
What  was  needed  was  a  device  that  would  give  reproduceable  readings 
for  a  given  suction  flow  rate  which  could  then  be  correlated  with  the 
appropriate,  approximate  tunnel  speeds.   Fine  adjustments  could  be 
obtained  in  other  ways  as  necessary.   At  this  point,  three  methods 
were  examined  for  determining  this  flow  rate:   an  orifice  meter,  a 
venturi  meter,  and  a  rotameter  or  flowrator.   The  first  method 
was  discarded  as  introducing  an  unwelcomed  large  pressure  drop 
(i.e.,  a  large  head  loss),  and  the  second  was  discarded  because  of 
the  limited  field  that  could  be  covered  with  sufficient  accuracy. 
The  third  method  might  have  been  rejected  on  the  same  grounds, 
except  that  for  the  rotameter,  various  size  floats  could  be  made 
to  cover  the  flow  rates  of  interest,  provided  the  bore  of  the 
instrument  was  sufficiently  large.   Most  importantly,  the  head  loss 
was  the  least  of  the  three  and  negligible  in  all  cases.   Although 
this  last  method  was  perhaps  the  most  costly  ordinarily,  a  Fischer 
&  Porter  flowrator  with  a  3"  bore  was  obtained  quite  inexpensively 
as  government  surplus  -  a  1948  NACA  item  now  considered  obsolete. 
This  instrument  was  traced  by  its  serial  number  to  obtain  the 
specification  sheet  which  indicated  that  it  had  been  built  and 
floats  designed  for  use  with  gasoline  and  kerosene.   Of  course  the 
calibration  curves,  even  if  they  were  available,  would  be  of  little 
use  since  it  was  now  intended  that  the  instrument  be  used  in  water. 
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As  the  end  result,  to  cover  the  upper  ranges  of  the  expected  flow 
rates  through  the  suction  system,  it  was  necessary  to  design  and 
build  one  float  in  addition  to  the  one  supplied  with  the  rota- 
meter.  Calibration  for  both  floats  was  performed  by  running 
water  through  the  rotameter  at  a  given  rate,  noting  the  meter 
reading  and  the  time  at  which  a  certain  weight  of  water  was  con- 
tained in  a  large  cylindrical  drum.   Note  that  this  weight  was 
greater  than  that  accumulated  in  the  time  for  reaching  a  steady 
flow.   The  water  was  allowed  to  continue  flowing  at  that  same 
rate  until  a  specified  greater  weight  was  in  the  drum,  and  the 
time  between  reaching  these  two  weights  was  recorded.   By  using 
the  elapsed  time  and  the  weight  of  water  accumulated,  the  flow  rate 
in  gallons  per  minute  was  obtained.   It  should  be  noted  that,  as 
expected,  these  data  plot  as  a  straight  line  for  both  floats. 
It  was  felt  that  although  the  results  obtained  represented  the 
flow  rate  to  within  a  few  gallons  per  minute,  the  curves  indicated 
the  more  important  result  that  the  readings  were  repeatable, 
the  primary  requirement  for  this  device.   As  long  as  a  given  flow 
rate  passing  through  the  rotameter  could  be  measured  consistently, 
then  for  a  given  approximate  tunnel  speed,  the  flow  through  the 
suction  system  could  be  set  near  the  proper  value  by  observing 
the  rotameter  reading. 
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With  a  rough  estimate  of  the  proper  amount  of  suction  obtained 
by  using  the  rotameter,  further  refinements  were  now  sought  which 
would  indicate  the  exact  suction  flow  rate  needed.   The  importance 
of  satisfying  this  condition  as  closely  as  possible  was  stipulated 
by  both  matching  with  the  theoretical  model  to  establish  the 
Blasius  type  boundary  layer  velocity  profile,  and  producing  as 
nearly  as  possible  the  zero  pressure  gradient  assumed  in  the 
analysis.   To  determine  the  exact  suction  needed,  three  methods 
appeared  as  likely  candidates:   the  injection  of  dye  at  the 
leading  edge  of  the  plate;  a  transverse  trace  made  with  hot  film 
anemometry  equipment;  and  finally,  the  distribution  of  pressures 
down  the  length  of  the  plate. 

The  first  method,  the  introduction  of  dye  into  the  flow  (an 
upstream  distance  of  1.25"  between  the  leading  edge  and  the  point 
of  injection  was  used)  entailed  injection  on  the  same  streamline 
that  would  "split"  the  leading  edge  if  the  flow  were  unaffected 
by  the  test  section  walls,  observing  visually  whether  it  was 
drawn  into  the  suction  channel  (too  much  suction)  or  passed  under 
the  plate  (too  little  suction).   Since  it  was  assumed  that  the 
flow  entered  the  test  section  parallelly,  this  streamline  was 
expected  to  lie  in  the  same  horizontal  plane  as  the  leading  edge; 
that  is,  the  injection  device  should  discharge  the  dye  at  the 
same  level  as  the  leading  edge.   The  device  for  accomplishing  this 
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was  required  to  be  placed  so  as  not  to  generate  turbulence  in  the 
free  stream.   It  was  believed  that  the  design  of  Mateyka  (9)  did 
not  take  this  factor  into  consideration.   Additionally,  it  was 
required  that  this  device  be  retractable  from  the  leading  edge 
when  the  suction  system  was  not  in  use  or  when  the  proper  amount 
of  suction  had  been  established.   Both  of  these  requirements  were 
satisfied  by  the  design  shown  in  figure  2.3.   Note  that  by 
inserting  the  device  from  the  top  of  the  nozzle  section,  any 
irregularities  caused  by  the  flow  past  the  shaft  would  be  drawn 
off  through  the  suction  system. 

In  the  second  method,  the  trace  made  by  the  hot  film 
anemometry  equipment  would  indicate  the  boundary  layer  velocity 
profile.   Correlation  of  this  profile  with  the  desired  Blasius 
type  profile  would  indicate  whether  the  proper  amount  of  suction 
was  being  used.   It  should  be  noted  that  this  was  never  expected 
to  match  exactly  due  to  the  inability  to  obtain  a  "perfect"  zero 
pressure  gradient. 

The  third  method  required  the  installation  of  pressure  taps 
down  the  length  of  the  plate  so  that  the  pressure  gradient  could 
be  monitored.  Figure  12  of  reference  (10)  indicated  the  erratic 
pressure  distribution  without  suction  and  compared  this  with  the 
distribution  attained  with  the  proper  amount  of  suction  at  the 
leading  edge  of  the  plate.   Based  on  these  results,  it  was  thought 
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Figure  2.3   Dye  injection  device 
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that  by  regulating  the  suction  so  as  to  obtain  the  proper  boundary 
layer,  the  pressure  distribution  should  appear  to  be  approximately 
constant.   To  monitor  these  pressures,  a  series  of  seven  pressure 
taps  were  installed  in  the  plate,  the  first  being  2-5/8"  from  the 
leading  edge,  with  each  succeeding  one  1-1/2"  downstream, 
alternately  staggered  1-1/2"  on  either  side  of  the  centerline. 
Care  was  taken  in  producing  the  holes  to  ensure  that  the  edges 
on  the  working  side  of  the  plate  were  made  as  squarely  as  possible 
by  drilling  from  the  reverse  side  and  then  lapping  the  plate.   At 
first  concern  developed  as  to  whether  the  effect  of  these  edges 
not  being  perfectly  square  would  introduce  any  error  into  the  system; 
for  certainly,  if  this  error  were  significant,  then  the  structural 
imperfections  from  one  hole  to  the  next  might  give  erroneous 
readings  and  explain  some  of  the  variation  noted  by  Mateyka  (9). 
Franklin  and  Wallace  (11)  in  evaluating  the  static  hole  error  for 
holes  about  the  size  that  were  needed  here,  explained  that  even 
with  a  hole  with  perfectly  square  edges,  a  residual  error  remains 
in  the  readings  due  to  the  local  change  in  the  boundary  conditions 
and  the  consequential  disturbance  of  the  boundary  layer:  and  so,  they 


obtained  a  plot  of  (error  in  pressure  reading  wall/shear  stress) 

^—      ,  where     Ur 


(C/T    )   versus   the   Reynolds  number    (Re   -    — K—      ,  where     11 


i 


is  the  friction  velocity  (Jr=  (fw/P  J1   »  d   is  the  diameter 
of  the  hole,  and  \)   is  the  fluid  kinematic  viscosity).   To  obtain 


23 

an  indication  of  the  size  of  the  error  for  this  problem,  the  shear 
stress  based  on  a  Blasius  velocity  profile  was  used  (6);  i.e., 


T  =  .332/uU, 


CO   N  ^OO 


a? 


v"  X 

Therefore,  for  a  free  stream  velocity  of  40  fps  at  60°F,  a  hole 
diameter  of  .040"  and  x  =  2-5/8",  the  distance  to  the  first  probe 
(these  being  the  conditions  that  would  introduce  the  greatest  error 
for  the  proposed  plate  design),  it  was  found  that, 


^p-  =    2IG.4  and  Tw=    8.44  x/O"3     Ib/^a 


From  the  (  C  /  Tw  )- Re  plot,  an  error  of  approximately  €=  .005904  psi 
or  .1634  in  H_0  was  found.   At  this  point  it  should  be  noted  that 
the  above  curve  was  based  on  data  obtained  with  holes  for  which 
unusual  care  (using  a  special  technique)  had  been  taken  to  ensure 
squareness  of  the  edges.   It  was  further  observed  by  Franklin  and 
Wallace  (11)  in  a  literature  survey  on  the  subject,  that  hole 
imperfections  resulting  from  the  method  of  hole  construction  for 
this  plate,  would  increase  the  error  up  to  as  much  as  five  times 
(i.e.,  up  to   e  ■  .0295  psi  or  .8171  in  FLO).   Upon  nondimensional- 
izing  this  latter  error  in  pressure  reading  with  the  dynamic  pressure, 
-L  ^  U^     ,  it  was  found  that  e/-L  o  U^    =  .002740.   By  com- 
parison, Mateyka  (9)  plotted  the  ratio  of  static  to  dynamic  pressure 
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and  reported  an  uncertainty  of  +  .005.   Note  that  the  error  found 
here  was  the  maximum  expected  for  any  flow  configuration.   It  was 
felt  that,  since  all  the  holes  were  identically  constructed  and 
only  the  relative  pressure  readings  on  the  plate  were  of  interest, 
the  static  hole  error  could  be  neglected. 

For  the  plate  designed  in  this  study,  a  comparison  is  made 
between  the  pressure  distributions  existing  with  and  without 
suction  (for,  as  an  example,   U  =  10.2  fps  and   7\  =  .701°)  is 

CO 

presented  graphically  in  figure  2.4,  and  is  discussed  more  com- 
pletely in  section  2.4C.   Basically  it  can  be  noted  that  variation 
of  the  suction  rate  does  not  significantly  alter  the  pressure 
distribution,  meaning  that  the  pressure  measurements  will  not 
provide  an  accurate  indication  of  the  exact  suction  needed  for 
a  given  free  stream  velocity. 

In  fact,  evaluating  the  methods  for  "fine  tuning"  the  suction 
rate  discussed  above,  the  following  order  of  utilization  evolves  - 
from  the  least  to  the  most  accurate: 

1.  pressure  measurements 

2.  dye  injection 

3.  boundary  layer  traverses  with 

anemometry  equipment 
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Note:   Successive  taps  alternate  1— "  either  side 
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Figure  2. A  Effect  of  leading  edge  suction  on  the  flat  plate 

pressure  distribution  for  ?\=   .701°  and  U^  =  10.2  fps. 
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2.4  Pressure  Gradient 

A.   Effect  on  flow  stability 

With  its  fixed  walls,  the  one  aspect  that  the  constant-cross- 
sectional-area  test  section  cannot  compensate  for  is  the  favorable 
streamwise  pressure  gradient  occurring  naturally  because  of 
boundary  layer  growth,  and  augmented  by  plate  insertion  into  the  flow 
(see  section  2.2A).   Therefore,  since  the  proposed  experimental 
configuration  assumes  a  zero  pressure  gradient,  before  implementing 
modifications  to  the  existing  facility,  it  is  essential  that  the 
effect  of  a  pressure  gradient  be  examined. 

The  strong  dependence  on  the  existing  pressure  gradient  of 
the  limits  of  stability  and  the  degree  to  which  a  small  disturbance 
would  be  amplified  have  been  shown  experimentally  by  Schubauer  and 
Skramstad  (7),  revealing  that  decreasing  downstream  pressures  cause 
a  damping  of  velocity  fluctuations  in  the  laminar  boundary  layer. 
Earlier,  Schlichting  (12)  theoretically  analyzed  for  stability, 
the  velocity  profiles  in  convergent  and  divergent  channels  resulting 
from  the  potential  flow  given  by  (J(*)  ■   Uoo  ~  (constant)  X 
with  a  similar  conclusion.   Later,  Schlichting  and  Ulrich  (13) 
carried  out  stability  calculations  using  a  modified  Pohlhausen 
approximate  method  and  determined  the  critical  Reynolds  number  as 
a  function  of  a  dimensionless  shape  factor,  defined  as: 

A  -  j£M 


27 

(where  for  the  purposes  of  the  ensuing  discussion,   o  is  the 
maximum  boundary  layer  thickness  on  the  plate  for  a  given  Uqu    , 
occurring,  of  course,  at  the  trailing  edge  of  the  plate).   This 
parameter  can  be  interpreted  physically  as  a  ratio  of  the  pressure 
to  viscous  forces. 

It  can  be  seen  from  a  plot  of  Schlichting  and  Ulrich  (13) 

i  i  c* 
of  the  critical  Reynolds  number   Re  +   =  — ^ —   versus  the 

%  6  c-nt     V 

shape  factor  j\    =  —  c±^l        ,  the  small  change  of  Yl   away 

^  cTx 
from  the  zero  pressure  gradient  case  for  this  favorable  one 

results  in  a  correspondingly  small  increase  in  the  critical 
Reynolds  number.   For  example,  if    _/\_   ■  .23,  an  apparently 
characteristic  value  by  subsequent  theoretical  computations  (see 
section  2.AB),  the  increase  in  the  critical  Reynolds  number  is 
approximately  90.   From  the  solution  of  Wazzan,  Okamura,  and 
Smith  (3),  such  an  increase  could  be  expected  from  only  a 
AT  =  Tw  -Tqo   -  3°F  (see  figure  3.1).   Since  it  is  expected 
that  experimental  observations  will  be  made  at  temperature  dif- 
ferences greater  than  this  (see  section  3.2A),  the  effects  of 
heating  the  plate  should  cause  a  more  significant  increase  of  the 
critical  Reynolds  number  than  the  existing  favorable  pressure 
gradient. 

B.   Design  estimate  of  the  magnitude  of  the  pressure  gradient 
For  the  particular  case  involved  here,  with  fixed  walls, 
boundary  layer  growth  on  the  sides  of  the  test  section  would 
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accelerate  the  free  stream  flow  as  if  it  were  a  perfect  fluid 
in  a  convergent  channel.   Since  the  thickness  of  the  layer  (and 
consequential  decrease  of  the  effective  free  stream  area)  is 
dependent  on  the  free  stream  velocity,  one  would  expect 
that  the  tunnel  speed  (in  addition  to  other  factors  such  as 
the  temperature  of  the  flow,  geometric  configuration  of  the 
test  section  -  it  does  naturally  converge  slightly  with  a 

difference  between  the  inlet  and  outlet  areas  of   Aj  -  A0   ■ 

2 
0.067  in  ,  etc.)  would  affect  the  magnitude  of  the  pressure 

gradient. 

Making  the  following  assumptions,  an  indication  of  the 
effect  of  the  pressure  gradient,  the  shape  factor  _/\_   ,  can  be 
calculated  approximately; 

1.   A  Blasius  laminar  boundary  layer  build-up  on  the  surface 
of  the  plate  which  implies  (6) 


SO,)-  5.o^ 

sVi  =  l-72f^f 


2.   A  fully  developed  turbulent  boundary  layer  on  the  three 
remaining  test  section  walls  (which  are  treated  as  flat  plates, 
i.e.,  corner  effects  occurring  due  to  finite  size  of  test  section 
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are  neglected)  where  the  boundary  layer  thickness  £   ,  based  on 
a  2%  velocity  defect,  increases  from  approximately  0.14  inches 
at  the  test  section  inlet  (for   U^  —  15  fps)  to  .25  inches  at 
the  outlet  (A).   It  is  further  assumed  that  the  values  measured 
at  the  nozzle  exit  (i.e.,  at  the  test  section  inlet)  are 
approximately  the  same  as  those  at  the  leading  edge  of  the  plate 
(the  point  on  the  plate  which  encounters  the  free  stream  first, 
and  from  which  all  measurements  relating  to  the  boundary  layer 
are  to  be  measured  in  a  streamwise  direction),  thereby  giving 
a  conservative  value  for  J\_ 

3.   Turbulent  boundary  layers  are  defined  by  the  1/7-th 

power  velocity  distribution  law  (25)  where 

i 


«.  =  (2Y- 


from  which  can  be  determined   5   =  o/Q     . 

4.  The  potential  velocity  gradient  is  linear. 

5.  The  plate  is  parallel  to  the  top  of  the  test  section 
and  inserted  to  a  distance  of  3/16". 

Since  the  displacement  thickness   0   ,  is  defined  as  the 
amount  by  which  the  external  potential  flow  is  displaced  by  the 
decrease  of  the  velocity  in  the  boundary  layer,  and  since  the 
potential  flow  distribution  essentially  determines  the  test 
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section  pressure  gradient,  by  calculating  the  difference  in 
"effective  cross-sectional  areas"  (caused  by  the  boundary  layer 
build-up  and  the  natural  convergence  of  the  test  section)  ,  the 
potential  velocity  gradient  can  be  determined.   Calculation  of  the 
shape  factor  follows  directly. 

For  example,  for  the  unheated  flat  plate  at  a  temperature 
of   "!«,-  60°F, 

40         .218 

Since  there  is  a  laminar  to  turbulent  boundary  layer  transi- 
tion at  the  nozzle  exit  between  the  tunnel  speeds  of  10  fps  and 
15  fps  (A) ,  it  is  difficult  to  ascertain  exactly  what  the  boundary 
layer  thickness  on  the  walls  of  the  test  section  will  be  at  the 
leading  edge  of  the  plate,  due  to  the  finite  distance  between  that 
point  and  the  nozzle  exit  (where  the  boundary  layer  thickness  is 
known).   Based  on  previous  assumptions  this  should  be  no  greater 
than  0.14  inches.   Assuming  that  this  value  may  be  used  for  a 
tunnel  speed  of  10  fps,  the  shape  factor  may  conservatively  be 
calculated. 

Um  (f  rQ  A 

10         .248 
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Notice  that  the  values  for  J\       are  quite  consistent,  under  the 

assumptions  that  have  been  made  here,  with  the  experimentally 

determined  maximum  value  of  Mateyka  (9),   _/Y  =  .233.   Note 

however  that  this  value  is  based  on  his  measurements,  conducted 

in  the  same  test  section,  where  the  highest  free  stream  velocity 

used  was  (J   ~  7.63  fps. 

C.   Experimental  measurements  and  procedures  for  reduction  of 
the  pressure  gradient 

Upon  initial  installation  of  the  plate  assembly  into  the  tunnel 

(i.e.,  at  zero  angle  of  incidence),  it  was  found  that  the  actual 

favorable  pressure  gradient  in  the  test  section  was  apparently 

significantly  larger  than  that  previously  estimated.   Possible 

explanations  were  advanced:   (1)  insufficient  suction,  but  it  is 

found  that  once  suction  is  applied  and  the  stagnation  region 

removed  from  the  leading  edge,  changes  in  the  suction  rate  do  not 

significantly  affect  the  pressure  gradient  (see  figure  2.4  for 

example);  (2)  misalignment  of  the  working  surface  of  the  plate 

relative  to  the  bottom  of  the  test  section,  but  measurements  later 

indicated  only  slight  convergence,  and  that  no  worse  than  the  area 

change  already  noted  in  the  preceding  section;  (3)  boundary 

layer  build-up  in  corners  which  was  not  adequately  accounted  for 

in  the  analytical  estimation  of  section  2.2B;  (4)  consideration 

of  the  pressure  reading  obtained  from  the  last  downstream  pressure 

tap,  yielding  misleading  results  for  the  pressure  distribution  over 

the  major  portion  of  the  plate. 
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The  significant  deviation  of  the  last  downstream  pressure  tap 
reading  (as  seen  in  figures  2.4,  2.5,  and  2.6)  is  attributed  to  the 
close  proximity  of  that  tap  to  the  trailing  edge  of  the  plate 
(.375").   It  is  felt  that  the  curvature  of  the  stream  as  it  seeks 
to  flow  around  the  square  edge  of  the  plate's  downstream  end  causes 
a  lowered  static  pressure  in  that  region.   Hence,  although  this  tap 
does  not  provide  data  consistent  with  the  essentially  linear  pressure 
gradients  obtained  from  the  remaining  locations,  the  data  so 
obtained  does  show  that  subsequent  adjustment  of  the  plate's  angle 
of  incidence  fails  to  correct  the  large  pressure  drop,  and  that 
small  pressure  gradient  data  must  be  taken  near  the  leading  edge 
or  central  portion  of  the  plate  so  as  to  remain  away  from  the  large 
drop.   Therefore,  comparing  experimentally  measured  with  the 
analytically/empirically  calculated  values  of  the  shape  factor, 
y\_    ,  presented  in  section  2.4B,  the  former  is  found  to  be 
larger  but  essentially  constant  for  all  free  stream  velocities. 
Without  considering  the  last  point,  on  an  average,   A    =  .8977 
whereas  consideration  of  the  last  tap  reading  gives  an  average 
value  of  J\      =   1.3616,  an  increase  of  approximately  52%. 

Of  the  four  explanations  advanced,  the  last  two  seem  to  be 
the  most  likely  for  the  disparity  between  the  calculated  and 
observed  pressure  gradient. 

At  this  point,  several  possibilities  were  considered  for 
reducing  the  test  section  streamwise  pressure  gradient: 
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Figure  2.5  Flat  plate  pressure  distribution  for  zero  angle  of 
incidence  to  the  mean  flow  at  various  flow  velo- 
cities. 
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Figure  2.6   Flat  plate  pressure  distribution  for  different 

angles  of  incidence  to  the  mean  flow  at  U^-   8.0  fps 
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(1)  modification  of  the  test  section.   Fix  within  the  test 
section  a  "flexible  false  bottom"  which  would  extend  from  within  the 
nozzle  section,  through  the  test  section  to  the  split  diffuser 
section,  the  contour  of  which  would  be  adjusted  by  deformation  at 
several  positions  along  the  length.   In  this  manner,  the  adjustable 
floor  contour  would  compensate  for  the  decreased  effective  area  for 
each  tunnel  speed. 

(2)  inclination  of  the  plate.   Utilizing  the  designed  versatility 
of  the  plate  assembly,  to  incline  the  working  surface  to  the  mean 
flow  without  drastically  altering  the  external  piping,  sealing  of  the 
plexiglass  frame  supporting  the  test  plate,  etc..  the  original 

1/32"  cork  gasket  sealing  between  the  plate  and  the  plexiglass 
frame  would  be  replaced  by  a  thicker,  pliant  rubber  gasket  material 
(3/8"  medium  density  closed-cell  neoprene  sponge  rubber  eventually 
was  chosen)  (see  figures  2.2  and  2.7).   In  this  way,  the  angle  of 
the  plate  could  be  adjusted  by  merely  compressing  the  gasket  more  at 
one  end  than  at  the  other  while  the  frame  and  piping  remained  fixed 
in  place. 

Either  method  of  correction  would  be  made  to  compensate  for 
boundary  layer  growth  on  all  four  sides  of  the  test  section  -  not 
just  that  on  the  plate  -  and  so  the  required  correction  would  be 
much  larger  than  that  of  only  the  displacement  thickness  of  the 
boundary  layer  on  the  plate. 
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Although  it  was  not  anticipated  that  the  second  correction  proce- 
dure would  completely  eliminate  the  pressure  gradient  as  perhaps  the 
first  might  have,  it  was  thought  that  the  reduction  would  be  signi- 
ficant enough  to  give  it  preference  over  the  more  mechanically  involved 
modification  of  the  test  section.   It  was  not  expected  that  the 
small  angle  of  incidence  would  significantly  affect  the  results. 

An  additional  advantage  claimed  by  the  second  procedure  can  be 
utilized;  it  would  require  that  the  plate  be  inserted  further  into 
the  free  stream  than  the  3/16"  calculated  in  section  2.2B.   Measure- 
ments of  the  turbulence  intensities  indicate  that  the  effect  of 
this  turbulence  is  discernable  further  into  the  free  stream  than 
the  measured  boundary  layer  thickness  (5)*  the  additional  distance 
between  the  leading  edge  and  the  top  of  the  test  section  gained 
by  the  thicker  gasket  would  ensure  that  the  "turbulent  intensity 
boundary  layer"  could  be  drawn  off  as  well  as  the  velocity  boundary 
layer.   This  point  will  be  clarified  in  subsequent  discussion. 

Since  the  experimental  analysis  herein  is  directed  toward 
confirming  that  the  designed  apparatus  functions  to  duplicate  a 
specified  model,  attention  is  focused  primarily  upon  establishing 
trends  -  that  is,  the  qualitative  rather  than  quantitative  behavior 
of  the  equipment.   Therefore,  no  excessive  effort  was  expended 
to  amass  a  huge  collection  of  data  for  different  suction  rates, 
angles  of  plate  inclination  to  the  flow,  etc.   Instead,  the  plate 
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assembly  was  operated  in  various  configurations  (see  for 
example,  figures  2.4,  2.5,  and  2.6)  to  see  how  closely  it  could 
satisfy  the  required  conditions,  then  assessing  the  significance 
of  any  disparity  between  the  two.   To  this  end,  the  following 
observations  are  made. 

There  are  two  factors  which  apparently  cause  the  appearance 
of  an  unfavorable  pressure  gradient  in  the  test  section:   (1)  no 
suction  at  the  plate's  leading  edge,  and  (2)  increased  inclina- 
tion of  the  plate  to  the  flow.   As  may  be  verified  by  figures  2.4 
and  2.6,  this  gradient  is  generally  small  with  respect  to  the 
total  head,  with  the  maximum  pressure  occurring  at  or  slightly 
upstream  of  the  plate's  midpoint. 

1.  It  appears  that  once  suction  is  provided  at  the  leading 
edge  (so  that  the  flow  no  longer  encounters  the  stagnation  region 
there  as  if  for  a  no-suction  configuration) ,  wide  variations  of  the 
suction  rate  affect  the  pressure  gradient  only  slightly  (see 
figure  2.4). 

2.  As  displayed  in  figure  2.6,  although  an  increased  plate 
angle  of  incidence  to  the  flow  does  tend  to  decrease  the  pressure 
gradient  initially,  further  inclination  after  a  certain  angle  will 
introduce  an  adverse  pressure  gradient  and  apparently  not  decrease 
the  magnitude  of  the  overall  gradient. 
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At  increased  free  stream  velocities,  the  effect  of  U^ 
on  the  pressure  gradient  seems  to  decrease;  that  is,  analysis  of 
figure  2.5  shows  that  the  value  of  \_  n  \\   z      decreases  at  a 

2  \    woo 

given  location  on  the  plate  as  U  increases. 

Also,  although  trends  in  the  pressure  distribution  were 
apparent,  appearing  consistently  in  the  limited  amount  of  data 
taken,  within  the  limits  of  the  accuracy  of  the  readings,  the 
pressure  gradients  for  a  free  stream  velocity  less  than  five 
feet  per  second  can  be  taken  to  be  zero. 

2.5  Design  Considerations  for  the  Plate 
A.  General  discussion 

The  details  of  design  for  the  plate  were  purposely  kept  as 
simple  as  possible  to  facilitate  construction  and  to  allow  versatility 
(see  figure  2.7  for  ensuing  discussion).   Guided  by  the  desire  that 
the  components  of  the  plate  assembly  be  replaceable  without  having 
to  rebuild  the  whole  assembly,  the  final  design  consisted  of 
three  parts:   the  actual  flat  plate,  a  plexiglass  supporting  frame, 
and  a  transition  section  which  would  channel  the  flow  from  a 
rectangular  cross-sectional  suction  channel  to  a  standard  size 
pipe  fitting. 

The  plate  would  fit  into  a  slot  in  the  frame,  would  be  held  in 
place  by  a  series  of  studs  embedded  in  the  plate,  and  would  be  sealed 
with  a  flat  gasket  "sandwiched"  between  the  two.   Pressure  tap 
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Piping  Transition  Section 


Plexiglass  Supporting 
Frame 


Sponge  Rubber 
Gasket 


Pressure  Taps 


S.S.  Studs 


Anodized  Aluminum 
Flat  Plate 


Figure  2.7   Exploded  view  of  the  parts  of  the  plate  assembly 
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fittings  also  embedded  in  the  plate  would  pass  through  holes  in  the 
frame.   All  such  holes  were  made  oversized  and  studs,  rather 
than  bolts,  were  used  in  securing  the  plate  to  the  frame  to  permit 
the  variable  longitudinal  inclination  of  the  plate  relative  to 
the  frame  required  for  the  pressure  gradient  elimination  discussed 
in  section  2. AC.   In  this  way,  plate  pieces  could  be  easily 
changed  without  disrupting  the  rest  of  the  apparatus,  sealing 
problems  would  be  minimized,  access  to  the  pressure  taps  would  be 
unhampered  by  a  top  plate  (i.e.,  there  would  be  no  problem  with 
having  to  seal  around  where  the  pressure  taps  come  out  of  the  plate 
due  to  the  gasket  seal).   Similarly,  for  the  magnet  chamber,  access 
to  the  magnets  would  be  unhampered  by  a  sealing  enclosure,  and  so 
virtually  any  reasonable  size  magnet  required  for  the  proposed 
electromagnetic  perturbing  mechanism  could  be  accommodated.   The 
dilemma  then  became  a  compromise  between  how  far  forward  the  magnet 
chamber  could  be  placed  to  the  leading  edge  while  still  allowing 
the  suction  channel  at  the  leading  edge,  and  the  requirement  that 
there  be  sufficient  structural  strength  between  the  two  so  that 
failure  would  not  occur  when  operating  at  maximum  tunnel  conditions. 

The  whole  apparatus  had  to  be  constructed  so  that  extraneous 
vibrations  of  the  plate  resulting  from  the  large  flow  rates  were 
kept  to  a  minimum.   To  this  end,  for  example,  rather  than  making  the 
leading  edge  of  the  plate  completely  planed  off  at  thp  angle  of 
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ten  degrees,  small  "lips"  were  left  on  both  sides.   Rather  than  the 
plate  fitting  flatly  on  the  frame,  the  latter  was  grooved  so  that 
the  plate  would  be  held  in  place  not  only  by  the  walls  of  the  test 
section,  but  by  the  recess  in  the  frame  as  well. 

In  designing  the  suction  channel,  a  number  of  factors  were 
considered.   To  utilize  the  testing  length  to  the  fullest,  the 
suction  channel  had  to  be  as  far  forward  as  possible,  thereby 
dictating  that  the  flow  be  drawn  off  vertically  (as  opposed  to  the 
more  favorable  small  angle  with  the  incoming  flow)  to  accommodate 
the  perturbing  mechanism  that  would  have  to  be  placed  on  the  plate. 
The  width  of  the  suction  channel  was  maintained  essentially 
constant  throughout  the  turn  to  ensure  that  the  flow  across  the 
inlet  remained  as  uniform  as  possible.   To  reduce  the  friction 
losses  with  the  relatively  large  flow  rates  that  would  be  experienced 

at  the  maximum  capacity  of  the  water  tunnel,  the  cross-sectional 

2 
area  of  the  suction  channel  was  increased  from  the  inlet  (A  =  1.69  in  ), 

2 
through  the  plate  and  transition  section,  to  the  piping  (A  =  4.26  in  ) 

that  carried  the  flow  to  the  pump.   As  the  flow  in  the  channel  would 

initially  be  turned  through  a  right  angle,  even  a  cursory  examination 

of  the  friction  losses  around  a  mitred  corner  indicated  the  advantage 

of  curving  the  channel  in  the  flow  direction. 
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B.   Selection  of  plate  material 

With  the  theoretical  design  considerations  taken  care  of,  there 
remained  the  problems  of  a  more  practical  nature  -  the  selection  of 
the  material  for  the  plate  and  then  the  design  of  the  plate  itself. 

For  purposes  of  this  investigation,  a  hard  anodized  6061  T6 
wrought  aluminum  plate  was  found  to  satisfy  all  the  requirements  for 
the  test  piece  (14  -  20).   Specifically,  since  no  large  localized 
stresses  were  anticipated,  a  softer,  more  easily  machineable 
core  material  could  be  used  under  a  hard  surface  coating.   The 
latter,  Al  O,  ,  while  strongly  adhesive  to  the  aluminum 
substrate,  provided  a  highly  dielectric,  corrosion  and  abrasion 
resistant  surface.   Both  materials  were  nonmagnetic  to  avoid 
shorting  the  magnetic  flux  path  needed  for  the  disturbance 
mechanism  as  discussed  in  chapter  4.   Further  motivation  for 
this  type  of  coating  stemmed  from  its  successful  application  on 
high  speed  airfoils  where  the  leading  edge  would  be  susceptible 
to  damage  due  to  the  surface  attrition  of  dust  and  rain. 

The  only  feature  of  the  anodic  coating  which  affected  the 
physical  specifications  for  the  plate  was  the  increase  in  each 
surface  dimension  by  about  one-half  the  thickness  of  the  anodic 
coating  produced.   Consideration  of  all  the  properties  of  the 
coating  was  felt  essential  to  ensure  that  the  plate  would  perform 
properly,  and  the  results  of  that  study  are  illuminated  in 
appendix  A. 
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C.   Surface  roughness 

Since  this  plate  was  designed  for  use  in  a  stability  study, 
the  effects  of  surface  roughness  and  waviness  on  transition  had  to 
be  considered  -  specifically,  the  fact  that  roughness  elements 
would  generate  disturbances  in  the  laminar  boundary  layer  in 
addition  to  those  already  present  there  due  to  the  flow  turbulence 
level.   If  the  former  disturbances  were  larger  than  the  latter, 
then  a  smaller  degree  of  amplification  would  be  required  then  for 
the  smooth  plate  surface  to  effect  the  transition.   Conversely,  if 
the  surface  were  very  smooth,  disturbances  generated  by  surface 
roughness  elements  would  be  less  than  those  already  present  due 
to  the  free  stream  turbulence,  and  the  presence  of  roughness 
would  have  no  effect  on  transition.   Therefore,  it  was  necessary 
to  determine  the  degree  of  surface  roughness  permissible  for  the 
flat  plate. 

In  the  limited  number  of  experiments  on  the  effect  of  distri- 
buted surface  roughness  on  transition,  the  results  have  been  some- 
what difficult  to  analyze  due  to  the  necessary  experimental  treat- 
ment of  the  leading  edge.   According  to  Lin  (21)  in  analyzing 
Holstien's  study  (22),  a  Reynolds  number  based  on  the  free  stream 
velocity,  the  nominal  roughness  height,  and  the  existing  viscosity, 
must  be  less  than  about  100  +  50  for  the  surface  to  be  considered 
aerodynamically  smooth.   The  experiments  conducted  later  by  Feindt 
(23)  indicated  that  the  grain  size,   K.  ,  in  the  presence  of 
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sand  roughness  had  a  significant  effect  on  the  critical  Reynolds 
number  only  after  the  value 


o  UcoK  -  12.0 


had  been  exceeded.   Prior  to  this  point  there  was  no  apparent 
change  in  the  critical  Reynolds  number,  but  after,  the  value  de- 
creased steeply.   For  this  problem,  using  maximum  operating  con- 
ditions,  K  >  334  microinches.   Hence  if  the  grain  size  were 
kept  less  than  this  value,  the  distributed  surface  roughness  should 
have  no  effect  on  transition.   Tests  on  the  surface  of  the  plate 
in  its  final  form  (including  the  surface  anodization  to  be  dis- 
cussed later)  were  conducted,  and  the  nominal  average  surface 
roughness  was  found  to  be  30  microinches,  certainly  well  within 
the  limits  discussed  above.   (This  would  correspond  to  Uco  ^s  —  8  1^ 
as  opposed  to  the  limit  of  120  noted  above.) 


CHAPTER  3 

THE  GENERATION  OF  HEAT  REQUIRED  TO  MAINTAIN  A  CONSTANT 
TEMPERATURE  FLAT  PLATE  IN  A  UNIFORM  PARALLEL  STREAM 

3 . 1  Selection  of  Heating  Method 

To  model  a  heated,  constant  temperature  flat  plate  at  zero 
incidence  to  a  parallel  laminar  flow,  a  special  method  for  heating 
this  plate  had  to  be  devised,  dependent  upon  assumptions  concern- 
ing the  temperature  dependence  of  flow  parameters  and  the  spatial 
dependence  of  thermal  properties  of  the  plate  (i.e.,  homogeneity 
of  material),  the  rate  of  surface  heat  flux  as  determined  by  the 
formation  of  the  thermal  boundary  layer,  the  transient  effects 
governing  the  time  required  to  attain  steady  state,  the  temperature 
difference  between  plate  surface  and  free  stream  and  associated 
buoyancy  effects,  and  so  forth.   Experience  and  theory  have  shown 
that,  due  to  the  large  temperature  gradients  near  the  leading 
edge  of  the  plate  caused  by  the  thinness  of  the  thermal  boundary 
layer,  the  rate  of  heat  transfer  for  a  given  temperature  would  be 
greatest  there,  decreasing  with  distance  down  the  plate  as  the 
boundary  layer  became  thicker.   Therefore,  regardless  of  any 
assumptions  made  in  trying  to  simulate  the  phenomenon,  if  a  con- 
stant temperature  were  to  be  maintained  on  the  surface  of  the  model, 
spatially  varying  amounts  of  heat  had  to  be  supplied  to  the  plate 
for  transfer  to  the  fluid,  with  the  greatest  amount  being  at  the 
leading  edge.   A  survey  of  the  literature  indicated  several  possible 
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methods  for  accomplishing  this  heating:   passing  a  vaporized  liquid 
(gas)  or  a  mixture  of  the  vapor  (gas)  and  liquid  on  the  reverse 
side  of  the  object  to  be  heated  or  cooled  (to  utilize  the  constant 
heat  of  condensation  at  the  appropriate  pressures)  either  over 
the  whole  surface  or  through  a  series  of  tubes  (much  like  a  heat 
exchanger)  (9,  10,  24);  the  use  of  a  number  of  rheostat  controlled 
strip  electric  heating  elements  (25,  26).   With  each  method,  there 
were  problems  in  establishing  and  maintaining  the  constant  tempera- 
ture.  For  example,  for  the  installation  under  consideration  where 
the  plate  would  be  mounted  horizontally  with  the  working  surface  on 
the  bottom,  heating  the  plate  by  the  condensation  of  steam  would 
cause  the  formation  of  film  or  drop-wise  condensation  on  the  upper 
surface  of  the  plate,  thereby  affecting  the  uniformity  of  the  heat 
transfer.  When  using  either  the  passage  of  a  fluid  through  a  series 
of  coils  or  the  electric  resistance  heating  of  a  series  of  heating 
elements,  the  uniform  heat  flux  on  the  working  surface  of  the  plate 
is  being  approximated  by  a  finite  distribution  of  sources,  so  the 
temperature  can  never  be  constant;  that  is,  there  will  always  be 
some  "waviness"  in  the  spatial  temperature  distribution  on  the 
plate.   However,  owing  to  the  success  of  several  experiments  (25, 
26)  where  heating  of  a  flat  plate  was  accomplished  by  use  of 
strip  heating  elements,  and  the  relative  simplicity  of  the  equip- 
ment needed,  this  method  was  chosen  for  this  experiment. 


3.2   Implimentatlon  of  Chosen  Method  -  Determination  of 
Distribution  of  Heating  Elements 

A.   Physical  limitations  affecting  design 

With  the  method  of  heating  now  chosen,  it  remained  to  actually 
design  the  plate  so  that  the  internal  heat  sources  would  closely 
approximate  the  actual  heat  transfer.   In  this  way,  the  spatially 
periodic  variation  of  temperature  on  the  working  side  of  the  plate 
caused  by  the  discrete  distribution  of  sources  separated  by 
adiabatic  regions  would  be  made  as  small  as  possible  within  the 
limits  of  the  physical  size  and  capacity  of  the  heaters.   For 
this  design,  a  temperature  difference,  A~T  ~  ~Kj      ~     lw   —   1°F 
along  the  length  of  the  plate  was  sought,  an  improvement  on  the 
"constant  temperature"  heated  flat  plates  used  previously  in  this 
water  tunnel  by  Mateyka  (9)  and  Morris  (10),  who  reported  a  varia- 
tion of  approximately  8°  -  15°C.   Note  that  in  formulating  this 
problem,  it  was  not  stipulated  that  both  the  surface  temperature 
and  heat  flux  be  specified  (the  two  conditions  that  should  match 
with  the  exact  laminar  thermal  boundary  layer  solution)  to  determine 
the  distribution  of  sources.   Since  the  equations  defining  conduc- 
tion in  the  plate  were  elliptic,  it  could  be  shown  that  the  imposi- 
tion of  Cauchy  conditions  on  the  surface  of  a  closed  domain  over- 
specified  the  problem,  and  for  an  open  domain,  made  the  solution 
unstable  (27). 
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Since  the  shape  and  size  of  the  plate  were  predetermined  by 
the  test  section  and  the  requirement  that  the  heated  plates  be 
interchangeable  (i.e.,  all  be  able  to  fit  into  the  frame  in  the 
test  section),  the  problem  was  reduced  to  determining  the  place- 
ment and  strengths  of  the  heaters.   Admittedly,  the  problem  might 
be  solved  with  very  little  computation  by  building  a  plate  with 
an  arbitrary  placement  of  the  heaters,  and  varying  the  heat  out- 
put of  each  until  the  surface  temperature  was  as  constant  as 
possible  for  that  distribution;  but  it  was  felt  that  an  analytical 
solution  prior  to  the  construction  of  the  plate  would  permit  the 
selection  of  spacing  and  distribution  of  strengths  of  the  sources 
to  better  simulate  a  constant  temperature,  and  the  "fine  tuning" 
could  still  be  performed  to  remove  any  minor  disparities. 

To  this  point,  it  has  been  assumed  that  the  amount  of  heat 
required  to  maintain  a  free-stream-to-wall  temperature  difference 
could  be  supplied,  but  no  indication  has  been  made  regarding 
its  magnitude.   Specifically,  the  maximum  temperature  difference 
that  must  be  maintained  is  required.   Since  the  tunnel  is  limited 
in  velocitv  and  the  plate  in  available  testing  length  (the  mini- 
mum being  the  position  of  the  perturbing  mechanism  as  measured  from 
the  leading  edge  of  the  plate)  to, 
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the  maximum  and  minimum  measureable  Reynolds  numbers   \\C€*       and 
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Figure  3.1   Plot  of  Re c  *   it   versus  wall  temperature  based  on  the 
results  of  Wazzan,  Okamura,  and  Smith  (3) 
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For  an  initial  estimate  of  the  temperatures  required?by 

plotting  Re  ^     versus  wall  temperature  from  the  data  of  Wazzan, 

a 

Okamura,  and  Smith  (2)  (figure  3.1),  and  then  superimposing  these 
limits  of  r\£c*  ,  it  can  be  seen  that  the  wall  temperature  must 

o 

remain  between  that  of  the  unheated  plate  and  that  with  a  difference 

'  w  ~oo   =  22.5°F.   For  design  purposes,  then,  a  maximum 

temperature  difference  of  25°F  is  used. 

B .   Analytical  solution  to  the  steady-state  conduction  equation 
in  a  rectangular  domain  with  rectangular  heat  sources 

Prior  to  a  mathematical  analysis  of  the  problem,  it  was  thought 
necessary  to  have  some  idea  of  the  design  in  mind.   The  one 
selected  for  this  case,  the  anodized  aluminum  plate,  is  the  same 
as  that  shown  in  figure  2.7  with  a  series  of  heating  elements 
embedded  in  it.   It  was  felt  that  the  degree  of  sophistication  in 
mathematically  modelling  the  plate  to  determine  the  distribution  of 
internal  heat  sources  required  for  a  given  heat  transfer  rate  from 
the  bottom  was  optional  to  a  certain  degree,  but  eccentric  geo- 
metries, assumptions  of  inhomogeneity  in  the  plate  in  either  the 
vertical  or  horizontal  direction,  and  the  like  each  only  made  the 
solution  more  complex,  but  not  necessarily  more  meaningful  for 
design  purposes.   For  purposes  of  this  analysis,  the  following 
assumptions  were  made: 


52 

1.  The  plate  is  of  rectangular  shape  with  length   31  , 
thickness  n  ,  and  width  w  ,  which  can  be  treated  two-dimens ion- 
ally  longitudinally  (  Bj»V\       and  w»  h  ). 

2.  The  heat  would  be  supplied  by  a  finite  series  of  heating 
elements  modelled  here  by  a  series  of  rectangular  sources. 

3.  The  top,  front,  and  back  of  the  plate  are  assumed  to  be 
insulated,  permitting  no  flow  of  heat,  and  the  bottom  with  the 
anodized  coating  acts  as  a  partial  insulator  -  a  resistance  to 
heat  flow.   The  efflux  of  heat  from  the  bottom  can  be  approximated 
by  the  Blasius-type  constant  wall  temperature  similarity  solution 
for  constant  property  boundary  layer  flow  (6); 


0  v   \>  Tx 


(3.2.1) 


Admittedly,  however,  the  assumption  of  constant  wall-to-free  stream 
temperature  difference  and  monotonically  decreasing  X  dependence 
implicit  in  this  formulation  would  not  be  compatible  with  the  tem- 
perature difference,  3      ,  obtained.   That  is,  due  to  the  discrete 
distribution  of  sources,  it  is  expected  that  some  oscillatory 
component  must  be  represented  in  the  solution,  the  amplitude  of 
which  would  decrease  directly  as  the  spacing  between  these  sources. 
Perhaps  the  heat  flux  would^  then  be  better  approximated  by 
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where 


represents  an  unknown  amplitude  of  oscillation  and  is  a  function 
of  the  spacing  between  successive  sources, 


and 

H(x-x-)  -  H(x-xj*0 

are  step  functions  restricting  a  particular  sinusoidal  contribution 
to  a  particular  interval.   Of  course,  the  exact  solution  would  require 
the  extensively  more  involved  solution  of  the  coupled  conduction  and 
convection  problem.   The  goal  of  this  analysis  is  to  reduce  the  amli- 
tude  of  these  oscillations  below  a  given  value,  and,  ideally,  the 
"corrective"  term  should  help  in  this  reduction.   For  purposes  of 
this  study,  however,  the  amplitude  is  assumed  to  be  small  (i.e., 

f ^   1),  so  only  the  leading  term  is  retained.   Modification  in 
strength  and  position  of  the  sources  is  made  until  this  assumption 
is  justifiable. 
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A.   The  plate  would  be  a  composite  consisting  of  aluminum 
of  thickness   D  and  thermal  conductivity  K^    ,    and  an  aluminum 
oxide  coating  of  thickness  C    and  thermal  conductivity  f\   , 
where  b »  C   and   K  »  K.   • 

5.  Each  region  is  homogeneous  (i.e.,  uniform  thermal  proper- 
ties) but  discontinuous  at  the  interface  between  the  two. 

6.  There  is  perfect  thermal  contact  at  the  interface. 
This  model  is  shown  in  figure  B.l. 

Before  proceeding  with  the  solution  to  the  problem  or  intro- 
ducing the  results  of  sample  calculations  for  various  source  con- 
figurations, a  few  observations  should  be  made  regarding  the 
validity  of  the  expected  solution.   As  noted  in  assumption  3  of 
this  section,  the  analysis  assumes  a  heat  flux  varying  inversely 
as  the  square  root  of  X  (equation  (3.2.1)).   Implicit  within  this 
assumption  are  (1)  a  constant  wall-to-free-stream  temperature 
difference  T^  - T oo  ,    and  (2)  the  requirement  that  there  be 

an  infinite  amount  of  heat  flux  as   X  *~^,  .   As  a  practical 

matter,  the  latter  statement  is  obviously  not  realistic.   In  fact, 
as  the  leading  edge  is  approached,  due  to  the  inability  to  physi- 
cally provide  sufficient  source  output  in  that  confined  space, 
the  wall  temperature  would  tend  toward  the  free  stream  temperature, 

T^j   *"  Tdo  •   That  is  to  say,  with  virtually  no  temperature 

difference  between  the  two,  the  heat  flux  would  tend  to  zero  in 


55 

that  region.   It  is  not  unexpected  that  this  mathematical  and 
physical  paradox  should  appear  when  keeping  in  mind  that  the  heat 
flux  used  for  this  model  is  derived  from  the  energy  boundary  layer 
equations  which  are  invalid  near  the  leading  edge  (   X  ^   0) .     \ 
Stating  this  conflict  in  another  way,  whereas  the  model 

requires  ~Liiv\        w    °°.         ,  physical  intuition  requires  that 

X  — "0+    f*~~ ' 

lion   lw  ~  * co  be  finite,  thus  violating  the  assumption  that  the  wall- 
to-free-stream  temperature  difference  be  constant.   In  short, 
the  mathematical  solution  in  the  region  near   X  =  0,  as  shown 
in  figures  3.2  and  3.3  for  example,  should  not  be  taken  as  being 
physically  meaningful. 

Additionally,  the  "hot"  regions  existing  in  the  neighborhood 
of  the  sources  and  the  relatively  "cold"  regions  between  them,  by 
having  a  larger  or  smaller  T^j —   l^  ,  respectively,  would  induce 
a  corresponding  increase  or  decrease  in  the  heat  transfer  in  those 
regions.   Unfortunately,  this  feature  cannot  be  incorporated  into 
the  present  analysis  without  the  simultaneous  solution  of  the  coupled 
problem,  as  pointed  out  in  assumption  3  of  this  section.   However, 
since  it  will  tend  to  decrease  the  amplitude  of  the  oscillations 
between  sources,  the  solution  obtained  herein  must  represent  the 
"worst  case"  (i.e.,  an  upper  bound  on  the  amplitude  of  the  oscilla- 
tion for  a  given  source  distribution). 
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Keeping  these  limitations  in  mind,  a  concise  restatement  of 
the  mathematical  model  is  made,  and  then  the  temperature  distri- 
bution is  sought.   Heat  is  generated  by  a  finite  series  of 
rectangular  heat  sources  in  one  region  of  a  two-region,  composite, 
rectangular  slab  which  is  insulated  on  three  sides,  but  which 
transfers  heat  by  convection  to  a  laminar  boundary  layer  on  the 
remaining  side. 

The  analytical  solution  is  obtained  by  first  determining 
the  proper  source  function  (or  since  the  temperature  gradients 
on  all  sides  are  the  prescribed  boundary  conditions,  the  Neumann 
function)  which,  when  found,  could  be  used  as  a  kernel  in  an 
integral  equation  incorporating  internal  sources  and  boundary 
conditions  to  define  the  temperature  distribution  within  the 
plate  (see  appendix  B,  sections  B.l  and  B.2  for  complete  deriva- 
tion).  Denoting  the  aluminum  as  region  "  31  "  and  the  anodized 
coating  as  region  "  C  »"  the  Neumann  functions,  N  ,  are  found 
to  be 

CO  CO  ,  I     — .    . 

TTrnv  />„e  IT^V 


*    *KB, 


(3.2.2) 
(3.2.3) 


'C 
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where  1)     tf0n  =  ±       ,   X  n  =  1   m  *  0 


(3.2.4) 


2)   Sources  are  assumed  to  be  in  region  1. 


3)   Eigenvalues   fb    are  determined  by 


Cot    ^nb      +  JV    CotfcnC    =  0 


(3.2.5) 


Assuming  that  the  heating  elements  could  be  represented  by  a  series 
of  step  functions 


s,0^)  =  I  SiR[H[x-(xt-4-)] 

H[x-(y;.A)]l|H[r(^4)]-H^ 

(  "^ 


(3.2.6) 


where 


.-tin 


■  strengths  of  the  l.    source 


^i  )Wt)   ■  position  of  the  center  of  the  I        source 
A   =  width  of  each  element 
P   =  thickness  of  each  element 
To  determine  the  eigenvalues  for  the  problem,  characteristic 
values  were  used  in  the  eigenvalue  equation  stated  above.   With 


b  =  .25 


C  =  .  oo2 


k  =  117  6<:u       «  -  S9  7  -^ 
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the  following  eigenvalues  were  found  for   (2>  ,    n  =  1,  2, 


&    =     3.1^830  1T       «      4tt 
.a    =     7.997673  ir        ~     8rr 


After  closer  scrutiny,  it  was  realized  that  £>n  was  approximately 

3n  ~  ^l?  n  ~    °»  1*  2'  * **  *    In  actuality*   the 

real  eigenvalues  were  less  than  this,  but  it  was  felt  that  the 


increasing  difference  between  the  real  and  approximate  values  was 
insignificant,  due  to  the  relatively  rapid  convergence  of  the 
solution  with  increasing   ft   • 

Therefore,  making  the  approximation  that   (3n  =  — i —  ,  the 
source  functions  became, 


^c  =  0  (3.2.8) 


where  now 
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This  significant  reduction  of  the  equations  implied  physically 
that  heat  transfer  through  the  anodized  coating  would  be  essentially 
one-dimensional  (normal  to  the  bottom  surface),  and  little  or  no 
enhancement  of  longitudinal  heat  transfer  could  be  expected  due  to 
the  presence  of  the  anodized  coating.   In  short,  the  plate  could 
be  treated  as  a  homogeneous  material. 

In  this  new  formulation  the  temperature  distribution  became 

fTnu-     r>     TTn  t 


_      2 


\ 


ab  K. 


s. 
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+  E 


cot 
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2b 
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ITinnX 
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•i.     Sin 
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m  =  1 
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(3.2.9) 
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a  a  2a  b  bJ 
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■rrri\2  1    /irm\/  rrn 


where 


2a.b      6a 


5,y^™b     ffi2ff«'f  d/ 


m=i 


(  /  ,  -  C  )  is  defined  by  equation  (3.2.1). 


C.   Asymptotic  series  solution  to  conduction  problem  as  an 
indication  of  the  required  source  distribution 

In  an  attempt  to  obtain  a  rough  indication  of  the  size  and 

position  of  the  sources  required,  one  can  examine  the  heat  flow  in 
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the  plate  on  a  still  simpler  basis  -  the  assumption  of  a  one-dimen- 
sional heat  flow.   The  justification  for  such  a  formulation  fol- 
lows . 

The  governing  equation  and  boundary  conditions  for  the 
rectangular  domain  are 


:i  +  a,  -  o 


where  (3.2.10) 

T.N-o)   =  Tw 

-&I(xa  0,a)=  0 


2T/  ,\       ^ 


|I(ro)^  ^fx)/Ki 


and    -c-  (w=b)  is  the  quantity  to  be  determined.   Q  (x)  is  again 
defined  by  equation  (3.2.1). 
Nondimensionalizing  with 

the  system  of  equations  (3.2.10)  becomes 
where 

©(To)=o         |§(ro)=ijM.        |®^i)  =  o 
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Further  assuming  that  0     can  be  expanded  in  the  parameter  (  b/  a  ), 


0=  a^tfe,  +(if©: 


-i-  .  • 


and  the  following  set  of  ordered  equations  evolve. 


^  a 

s   2. 

^/2 


-  0 


-f 


-  O 


and  so  on. 

Solving  for  each  order  and  applying  boundary  condtions  away  from 
the  end  points,  it  is  found  that  to  second  order, 


^4 


+  1- 

12.8 


-^"-i^ 


+  .  •• 


(3.2.11) 


Notice  that  this  solution  is  valid  as  long  as  no  attempt  is 
made  to  satisfy  boundary  conditions  at  either  end.   This  suggests 
the  presence  of  a  boundary-layer  type  solution  or  a  singular  per- 
turbation problem,  an  observation  substantiated  by  Van  Dyke's  rule 
for  the  uniformity  of  a  solution  in  a  given  domain  (28) ;  "A  per- 
turbation solution  is  uniformly  valid  in  the  space  and  time 
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coordinates  unless  the  perturbation  quantity  is  the  ratio  of  two 
lengths."  Note  further  that  in  solving  the  problem  thus  far, 
the  basic  equations  have  been  changed  in  nature  from  elliptic  to 
parabolic,  thereby  permitting  the  imposition  of  Cauchy  boundary 
conditions  on  one  boundary  (see  section  3.2A)  to  see  what  the 
conditions  should  be  on  the  opposite  one. 
From  equation  (3.2.1) 

b@  _     lo  9  f  x) 


**)         KJT^-To) 


st 


1  "8  J 


(3.2.12) 


+ 


-41 


4- 


35 


so 


8  J7-        /28 


(3.2.13) 


rj    =  1,   and 


Evaluating  the  solution  at  the  top  of  the  plate, 
realizing  that  except  very  close  to  the  leading  edge, 
a  good  representation  of  the  required  heat  flux  needed  at  the 
top,   Q(  X  ,  b  ),  is  given  simply  by, 

Q(\b)--  £<x) 

since    f  k-\  «   1 
To  represent  the  continuous  heat  flux  distribution  by  a  discrete 
distribution  of  heating  elements,   Q  (  x  »  b  )  is  integrated  over 
some  range,  AX  ,  which  is  determined  by  the  size  and  output  capa- 
city of  a  heating  element.   That  is 
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l(4*,b)=  .U4-  Ka  ffr  {%-  (T, - T«)(f  xc  +  ax1  -  f*T) 


(3.2.14) 


D.   Discussion  of  plate  temperatures  for  different  source 
distributions 

To  assess  the  capability  of  various  distributions  of  heat 
sources  for  providing  a  required  temperature  distribution  on  the 
plate,  two  examples  are  presented  as  sort  of  "numerical 
experiments."  For  both,  it  is  assumed  that  twelve  heaters  supply 
the  required  heat  output,  and  further,  that  each  can  individually 
provide  any  specified  amount.   Additionally,  although  three  variables 
x^ould  ordinarily  be  available  for  adjustment  (i.e.,  horizontal 
and  vertical  spacing,  and  source  strength)  due  to  the  limitation  on 
the  vertical  spacing  by  the  thinness  of  the  plate,  only  the  heater 
strength  and  horizontal  spacing  are  assumed  adjustable.   Physical 
parameters  applicable  to  both  examples  are  presented  in  table  3.1. 

The  first  example  specifies  a  distribution  of  constant 
strength  heaters  with  variable  spacing  (such  as  used  in  reference 
(26)).   This  spacing  is  determined  by  finding  the  total  heat  flux 
required  per  heater,  and  then  using  the  one-dimensional-analysis, 
equation  (3.2.14),  to  estimate  the  distance  along  the  plate  re- 
quired to  transfer  that  much  heat  by  convection.   The  heater  is 
then  assumed  to  be  placed  in  the  middle  of  this  interval.   Strengths 
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3  =  12.25" 

b  =  .25" 

C  =  .002" 

K  =  .340  Btu/hr   ft    °F 

K  =  117  Btu/hr   ft    °F 

Kc  =  .597   Btu/hr   ft    °F 

?r  =  8.03 

Uco  =20.0  fps 

AT  =  25.0°F 

n  =  .001 

\)  =  1.22  x  10"5   ft2/sec 

A  =  .25" 

E  -  .015" 

VS.  -  95,197.00818  l 
Ij>    La 


u 


ta  in3 


H  =  50.9983 


watts 

m3/2 


Table   3.1 


65 

and  positions  of  the  heat  sources  for  this  example  are  listed  in 
table  3.2.   Note  that  the  first  three  heaters  are  not  spaced  in 
accordance  with  equation  (3.2.14)  due  to  their  width.   That  is, 
the  theory  would  require  that  they  be  spaced  successively  at 

X  =  .0425",  .2125",  and  .5525".   With  the  heater  width  speci- 
fied as  Q   =  .25",  such  spacing  is  impossible,  so  they  are 
positioned  as  closely  as  physical  dimensions  will  permit.   The 
resulting  temperature  distribution  is  plotted  as  a  percentage 
of  the  wall  to  free  stream  temperature  difference  in  figure  3.2. 

The  second  example  specifies  a  distribution  of  variable 
strength  heaters  with  constant  spacing  (such  as  used  in  reference 
(25)).   The  required  strengths  are  found  in  much  the  same  manner 
as  for  example  1,  except  that  here,  a  distance  along  the  plate 
is  specified  and  the  heat  transferred  by  convection  in  this 
interval  evaluated.   The  appropriate  source,  in  turn,  is  assumed 
to  supply  this  quantity  of  heat  and  to  be  positioned  in  the  center 
of  the  interval.   Table  3.3  lists  the  source  strengths  and  posi- 
tions for  this  example.   Temperature  variations  over  the  entire 
length  of  the  plate  as  well  as  details  of  this  variation  in  a 
characteristic  interval  are  seen  in  figures  3.3  and  3.4, 
respectively. 

Discussing  these  two  examples  qualitatively,  with  reference 
to  figures  3.2  and  3.3,  it  is  noted  that  both  formulations  reveal 
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Source 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 


Strength  (S-A)    P^) 

in 


7,933.0840 
7,933.0840 
7,933.0840 
7,733.0840 
7,933.0840 
7,933.0840 
7,933.0840 
7,933.0840 
7,933.0840 
7,933.0840 
7,933.0840 
7,933.0840 


X.(ln) 

Ljfcttn) 

.1250 

.2 

.3750 

.2 

.6250 

.2 

1.0630 

.2 

1.7440 

.2 

2.5950 

.2 

3.6150 

.2 

4.8060 

.2 

6.1680 

.2 

7.6990 

.2 

9.4000 

.2 

11.2720 

.2 

Table  3.2 
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Source        Strength  (^-a  )  ( ^-) 


in- 


1  27,199.0980 

2  11,266.2350 

3  8,644.8866 

4  7,287.9764 

5  6,420.8362 

6  5,804.8797 

7  5,338.1377 

8  4,968.6176 

9  4,666.6274 

10  4,413.8061 

11  4,198.1029 

12  4,987.6405 


XjCln) 

#ln) 

.5 

.2 

1.5 

.2 

2.5 

.2 

3.5 

.2 

4.5 

.2 

5.5 

.2 

6.5 

.2 

7.5 

.2 

8.5 

.2 

9.5 

.2 

10.5 

.2 

11.5 

9 

Table  3.3 
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Figure  3.2   Plate  surface  temperature  distribution  for  constant 
strength  sources  with  variable  spacing 
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Figure  3.3   Plate  surface  temperature  distribution  for  variable 
strength  sources  with  constant  spacing 
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Figure  3.4   Detail  ofthe  flat  plate  surface  temperature  distribution 

in  a  characteristic  interval  for  the  example  of  a  variable 
strength  -  constant  spacing  source  distribution 
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inconsistencies  in  the  vicinity  of  the  leading  edge  (as  expected 
from  the  discussion  in  section  3.2.B),  higher  temperatures  in 
regions  of  the  heat  sources  with  lower  temperatures  in  the  adia- 
batic  regions  between  successive  sources,  and  larger  amplitude 
temperature  variations  for  greater  heater  spacing. 

On  the  other  hand,  whereas  in  example  1,  the  *-emperature  fluc- 
tuations increase  in  a  streamwise  direction  due  to  the  increased 
lateral  heater  spacing,  in  example  2,  they  decrease  because  of 
the  reduced  heat  transfer.   Additionally,  the  average  plate 
temperature  appears  to  constantly  increase  for  example  1,  but 
remains  essentially  constant  over  a  major  portion  of  the  plate 
for  example  2.   Unlike  the  results  of  the  former,  where 
the  temperature  variation  is  fairly  symmetrical  around  successive 
sources,  the  ununiform  spatial  temperature  oscillations  for  the 
latter  are  thought  to  be  caused  by  lateral  heat  transfer  in  the 
plate  arising  because  of  the  variable  source  strengths.   Note  from 
the  amplitude  irregularities  shown  in  figure  3.4  that  the 
horizontal  spacing  between  adjacent  heaters  is  apparently  not 
the  only  criterion  for  determining  the  amplitude  of  the  temperature 
variations. 


CHAPTER  A 

GENERATION  OF  SMALL  AMPLITUDE  WAVE  DISTURBANCES 
IN  THE  FLAT  PLATE  BOUNDARY  LAYER 

4 . 1  Selection  of  Generation  Method  Used 

To  select  the  means  for  generating  oscillatory  disturbances 
in  the  boundary  layer  of  the  designed  flat  plate  assembly,  it  is 
specified  that  the  disturbance  mechanism  be  capable  of  producing 
both  "finite"  amplitude  disturbances  and  the  "infinitesimal 
Tollmlen~Schlichting  waves"  assumed  in  conventional  stability 
analysis.   This  latter  capability  should  permit  a  reasonable  com- 
parison between  experimental  results  and  an  analytical  solution 
to  the  problem  under  consideration  (references  (2)  and  (3)). 

A  number  of  different  methods  have  been  tried  in  similar 
experiments.   Schubauer  and  Skramstad  (7),  concerned  with  sub- 
sonic speeds  and  frequencies  of  10  -  260  Hz  ,  investigated  the 
effects  of  sound,  first  by  installing  a  loudspeaker  at  the 
leading  edge,  and  then  by  introducing  the  sound  into  the  boundary 
layer  through  a  small  hole  in  the  flat  plate;  but  they  found  that 
resonance  effects  and  the  complexity  of  the  resulting  wave  patterns 
excessively  complicated  boundary  layer  perturbation  with  controlled 
sound  intensity.   They  further  used  successfully  an  electro- 
magnetically  vibrated  ribbon  stretched  parallel  to  the  surface  of 
the  plate  and  to  the  leading  edge.   This  method  has  since  been  used 
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by  other  authors  ((29)  and  (30),  for  example).   Laufer  and  Vrebalo- 
vich  (31),  concerned  with  the  laminar  supersonic  boundary  layers 
and  frequencies  of  5000  -  50,000  Hz  »  also  tried  using  sound  waves, 
in  this  case  generated  by  a  cutout  in  a  second  flat  plate 
(oriented  parallel  to  the  main  flat  plate);  but  the  radiated 
field  proved  to  be  too  wide,  disturbing  the  boundary  layer  along 
the  whole  length  of  the  plate.   They  tried  an  electromechanical 
device  where  a  thin  diaphram  mounted  flush  with  the  surface  of 
the  flat  plate  would  be  made  to  oscillate  by  acting  as  one  plate 
of  a  capacitor;  but  they  found  that  insulation  break-down  in  the 
capacitor  at  low  voltages  and  shielding  of  their  hot  wire  amplifiers 
from  the  high  frequency  voltage  source  radiation  field  made  this 
method  untenable  as  well.   A  third  (successful)  method  used  by  them 
was  essentially  a  high  speed  valve,  where  a  narrow  slit  in  the 
plate  would  be  periodically  opened  and  closed  to  permit  an  air  pulse 
of  any  required  strength  and  frequency  to  be  injected  into  the 
boundary  layer.   Polymeropoulos  and  Gebhart  (32)  and  later 
Dring  and  Gebhart  (33)  perturbed  a  natural  convection  boundary 
layer  with  a  wire  ribbon,  but  they  oscillated  the  frame  on  which 
the  ribbon  was  stretched  parallel  to  the  face  of  the  plate  rather 
than  the  ribbon  directly. 

After  reviewing  all  of  the  methods  used  above,  the  vibrating 
ribbon  was  chosen  for  this  experiment.  In  essence  the  sinusoidal 
oscillation  is  obtained  by  passing  an  alternating  current  of  a 


74 


desired  frequency  through  the  ribbon  in  the  presence  of  a 
constant  magnetic  field.   For  the  forces  to  be  normal  to  the 
plate,  the  magnetic  field  has  to  be  perpendicular  to  and  in 
the  same  plane  as  the  undisturbed  wire  ribbon.   To  produce 
such  a  field  ,  a  bipolar  or  horseshoe  type  magnet  is  used  with 
the  wire  located  equidistant  from  the  two  poles.   An  alternate 
approach  might  be  to  run  a  constant  current  in  the  wire  and 
oscillate  the  magnetic  field;  but  it  was  felt  that  this  large 
alternating  magnetic  field  would  most  probably  induce 
voltage  fluctuations  in  the  probes  used  to  detect  the  stability 
of  the  flow.   Therefore,  the  first  method  was  used. 

4 . 2  Compatibility  of  Experimentally  Generated  Wave 
with  the  Assumed  Tollmein-Schlichting  Wave 

At  this  point,  the  question  might  be  raised  as  to  whether 
such  a  device  could  closely  enough  model  the  two-dimensional 
disturbances  assumed  in  conventional  stability  analysis.   Gaster 
(34)  mathematically  shows  that  "a  ribbon-type  disturbance  ex- 
cites a  travelling  wave  system  downstream  of  the  source.   This 
mode  will  have  a  wave  number,  o^   ,  identical  to  that  of  the  free 
eigenmode  which  has  the  same  frequency  as  the  ribbon."  The  complex 
nature  of  this  wave  number  confirms  that  "the  mode  excited  by 
the  vibrating  ribbon  is  indeed  one  having  spatial  growth  as 
observed  in  experiment." 
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Lin  (35),  in  analyzing  the  work  of  Schubauer  and  Skramstad  (7), 
points  out  the  conformity  of  the  free  oscillations  observed  by  the 
latter  with  the  two-dimensional  theory,  noting  that  these  oscilla- 
tions are  known  to  be  three-dimensional  and  should  be  expected  to 
be  such  due  to  their  instability  at  sufficiently  high  Reynolds  num- 
bers.  He  shows  that  the  two-  and  three-dimensional  disturbances 
have  nearly  the  same  second  branch  point  for  the  neutral  stability 
curve  and  so  are  indistinguishable,  concluding  that  the  oscillations 
in  the  boundary  layer  become  progressively  more  three-dimensional 
as  the  Reynolds  number  of  the  flow  increases.   Similar  conclusions 
have  been  drawn  by  Klebanoff  and  Tidstrom  (30).   So,  it  would  seem 
that  any  attempt  to  generate  only  two-dimensional  disturbances  (as 
assumed  for  the  theoretical  model  in  this  analysis)  is  impossible, 
but  this  difficulty  had  been  removed  earlier  (at  least  for  an  un- 
heated  plate)  by  Squire  (36)  who  proved  that  a  two-dimensional, 
incompressible  flow  becomes  unstable  at  larger  Reynolds  numbers 
when  a  disturbance  within  it  is  assumed  to  be  three-  rather  than 
two-dimensional.   Hence  the  critical  Reynolds  number  is  determined 
by  the  latter  type  of  disturbance. 

A. 3  Determination  of  Significant  Parameters  Affecting  Vibrating 
Ribbon 

A .   Explanation  for  analytically  modelling  the,  vibrating  ribbon 

For  design  purposes,  it  is  thought  necessary  to  obtain  some 

sort  of  analytical  solution,  prior  to  the  construction  of  the 

perturbing  mechanism,  to  resolve  two  problems.   First,  since  such 

a  device  has  been  used  successfully  for  forced  convection  studies 
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exclusively  in  air,  where  one  would  expect  the  viscous  damping  to 

be  much  smaller  than  that  in  water,  prior  to  use  in  the  latter 
medium  an  indication  of  the  effect  of  viscous  drag  on  the  smooth 
oscillations  of  the  ribbon  should  be  examined  -  specifically,  the 
phase  and  form  distortion  of  the  generated  wave  and  change  in  its 
amplitude.   Nondimensionalization  of  the  governing  equation  with 
respect  to  significant  parameters  dispels  the  hope  that  the  relative 
smallness  of  the  coefficient  of  the  damping  term  would  permit  this 
term  to  be  completely  excluded  in  an  approximate  solution.   There- 
fore, the  nonlinear  viscous  damping  term  is  retained  in  the  formula- 
tion of  the  problem.   Secondly,  it  is  desired  to  obtain  an  idea  of 
the  magnitude  of  the  driving  force  needed  to  vibrate  this  ribbon. 
For  these  reasons  the  following  analytical  problem  is  posed. 

B.  Analytical  solution  of  the  nonlinear,  fourth  order  equation 
defining  the  motion  of  a  "beam"  string  with  viscous  damping 
being  driven  by  an  electromagnetic  force 

To  identify  all  the  physical  processes  taking  place  and  to 

facilitate  the  formulation  of  the  mathematical  model,  a  statement 

of  the  problem  to  be  solved  is  made.   A  wire  ribbon  of  length  1, 

width  Cr   ,  and  thickness  K  ,  is  stretched  parallel  to  a  flat  plate 

at  a  close  distance  and  constrained  at  each  end  with  built-in  end 

conditions  (that  is,  slope  and  deflection  at  end  points  are  zero). 

It  is  forced  to  vibrate  electromagnetically  in  a  viscous  medium 

by  passing  an  alternating  current  through  the  ribbon  of  frequency 

(-3    in  the  presence  of  a  constant  magnetic  field. 
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In  deriving  the  equation  for  the  motion  of  the  disturbance 
ribbon,  the  following  assumptions  are  made: 

1.  Motion  takes  place  in  only  one  plane,  and  in  this  plane, 
each  particle  moves  in  purely  a  transverse  direction  (i.e.  ,  at  right 
angles  to  the  equilibrium  position  of  the  ribbon). 

2.  Deflection  of  the  ribbon  is  so  small  at  all  times  that 
its  change  in  length  when  displaced  has  no  effect  on  the  tension. 
That  is,  the  slight  rotation  of  the  cross-sections  and  elongations 
of  the  ribbon  are  negligible  so  the  effects  of  structural  damping 
(37),  (38)  and  the  spring  restoring  force  are  neglected. 

3.  Since  the  maximum  amplitude  of  the  ribbon  is  assumed  to 
be  small,  for  simplicity  of  the  solution,  it  is  further  assumed 
that  there  is  no  change  in  the  magnetic  field  over  the  range  of 
displacement. 

4.  Viscous  drag  coefficient  of  the  ribbon  is  essentially 
constant,  representing  almost  exclusively  form  drag.   Because 
the  length  of  the  ribbon  is  much  greater  than  the  width,  which 

in  turn  is  much  greater  than  the  thickness,  it  can  be  assumed  that 
this  ribbon  acts  as  an  infinitely  long  flat  plate  moving  normal 
to  a  given  flow.   This  coefficient  has  been  determined  to  be  a 
constant,  independent  of  the  speed  of  the  flow  (39): 

CD  -  2.0 
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5.  The  ribbon  is  treated  as  a  "beam"  string  to  account  for 
the  bending  forces. 

6.  Spatial  boundary  conditions  require  that  both  ends  be 
clamped  in  place:  that  is,  the  slope  and  displacement  are  assumed 
to  be  zero  for  all  time.   To  establish  temporal  boundary  condi- 
tions, it  is  stipulated  that  the  effect  of  damping  will  be  a 
maximum  when  the  velocity  of  the  ribbon  is  the  greatest  and  zero 
when  the  ribbon  is  at  rest. 

With  these  restrictions  in  mind,  the  equation  describing  the 
forced  oscillation  is  derived. 
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with  spatial  boundary  conditions 


l£'=f=0       at      F 


O 


Temporal  boundary  conditions  are  to  be  established  shortly. 
Initial  significance  attached  to  multiple-time  scaling 
was  soon  revealed  as  unfounded,  and  so  the  solution  to  equation 
(A.  3. 3)  was  obtained  by  a  double  expansion  in   (2>  (associated  with 
the  viscous  damping)  and   £  (delineating  the  relative  importance 
of  bending  rigidity  and  ribbon  tension).   (See  appendix  C  for 
further  elucidation).   Hence  the  nondimensional  amplitude  can  be 
written  as 

I  =  [  g*(f>r>^MfA)  +  •••] 

(4.3.4) 

+ 


»    •   • 


where  essentially  the  functions  p1  and  <  represent  the  first 
order  corrections  to  the  forced  oscillation  of  a  massless  string 
due  to  the  effects  of  viscous  damping  and  bending,  respectively. 
Substituting  the  assumed  form  of  j   into  equation  (4.3.3)  an 
ordered  set  of  equations  in  powers  of  p      results  for  zeroth  order 
in  €  .   Imposing  the  spatial  boundary  conditions  and  applying 
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the  intricate  matching  conditions  in  regions  near  the  ribbon's 
end  points  (where  it  is  confirmed  that  the  uncorrected  solution 
must  satisfy  only  the  homogeneous  spatial  boundary  conditions  - 
that  is,  temporal  boundary  conditions  need  not  be  considered 
until  the  effects  of  damping  are  introduced),  the  zeroth  order 
solution  in  both  ft  and  £  is  found  to  be 


(f/r)  =    Xz   tan  |  sin°^  +  cos<*J  -1 


5LO  T       (A. 3.5) 


Note  again  the  Pq    is  expected  to  be  the  major  contribution  to 
the  solution  for  the  ribbon's  displacement. 

With  po     thus  determined,  the  constant  coefficients  in 
the  expression  for  p  1  (  F  ,  T  )  can  be  found  by  satisfying 
the  matching  conditions  and  by  realizing  that  the  viscous 
damping  must  be  zero  at  the  time  when  the  local  velocity  of  the 
ribbon  is  zero.   Therefore, 


+  K/f)  cos5t  +  K  Cos7r  +... 
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(4.3.6) 


where 
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K,=  irS  (tan^af^f   t-    cos  «f  -  lj 
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J£       /  "fc^n- ^  Sin  3«(  *c   4-  Cos  3*  £  -  1 


t/      =  K,    /tan^z  sinS^-h  cos  5*^  -1 

^3  875  V  tannin  ^    i-   ^s^^  -  1 


K 


5'45  ^  ta.n-^-5i"noi|    +-  Cos  <^  -  1 

Again  applying  matching  conditions,  it  can  be  shown  that  for 
first  order  matching  in  £   ,   "0  (£)  =  ]  £  and  so 


K0(|)t)=  ^Uo|  -tany  sm^  -  coscx^ 


SinT 


(4.3.7) 


Extending  this  procedure  further,   l(   and  additional  correction 
terms  can  be  obtained,  but,  for  purposes  of  this  analysis  where 
even  the  first  order  corrections  are  shown  to  be  small,  it  is 
felt  that  higher  order  approximations  are  unnecessary. 

In  an  attempt,  now,  to  obtain  some  idea  of  the  effect  of 
the  viscous  damping  and  bending  on  the  ribbon's  movement, 
typical  values  for  the  problem  are  used  to  evaluate  the  solution 
(see  table  4.1). 
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Material:   spring  phosphor  bronze  ribbon  (94.75%  Cv    ,  5%  5n  ,  .25%  P) 

f£  =10  lb. 

K  =  .0035  in 

K  -  .006  in 

0*  =  .0625  in 

m  =  .37  x  10~A  slugs/ft 

Pr  =  17.05  slug/ft3 

F  =  16  x  106  psi 

1  -  9  in 


K3 

1 

C     =  2.0 


I     ■  — r-r—  =  6.51  x  10    in 
o        12 


^o 


=  1.9A  slugs/ft3 


Table  4.1 
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The  significant  parameters  of  the  problem  can  be  determined: 

6  -  2.0(o  x  /0'5 

0^  =  1.765 

6  =  .  23<?  ■ 

X  is  postulated  to  be  of  unit  order. 
Hence,  the  time  dependent  amplitude  of  the  ribbon  at  its  midpoint 
is 

■f  =   [  ,  \530  Sin  r     -  {.00073    Cost    -.000  21  cos3l 

+  .0000]   C05  5r  +-. .  A  + ..  ♦  1 
-     .00340  smr  +    (?(^f£)] 

From  this  result,  it  is  seen  that  the  effect  of  the  damping 
is  to  cause  a  damped  ribbon  to  lag  an  undamped  one  (i.e.,   6=0 
solution).   Additionally,  since  the  ribbon  is  required  to  complete 
one  cycle  of  its  oscillation  in  the  same  period  as  that  of  the 
driving  force,  it  is  verified  that  the  damping  will  tend  to  speed 
up  the  ribbon  In  some  regions  (immediately  after  the  maximum 
displacement  is  reached)  forcing  it  to  go  slower  than  the  undamped 
one  in  others  (near  regions  of  zero  displacement).   These  two 
expected  effects  are  qualitatively  displayed  in  figures  4.1  and 
4.2;  quantitatively,  the  effects  are  much  smaller  than  those  shown 
as  substantiated  by  the  above  sample  calculation. 
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Figure  4.1   Effect  of  damping  and  bending  on  the  displacement 
of  a  vibrating  ribbon. 
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Figure  4.2   Effect  of  damping  and  bending  on  the  velocity  of 
a  vibrating  ribbon. 
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C.   Simplified  solution  to  the  same  problem 

Since  it  has  been  shown  in  the  preceding  section  that  the 
effect  on  the  wave  form  of  a  vibrating  ribbon  due  to  damping  and 
bending  is  small,  for  subsequent  analysis,  higher  order  correc- 
tions are  neglected  and  the  displacement  is  approximated  by 
the  zeroth  order  solution.   That  is, 


(4.3.5) 


sm  r 


or  expressing  all  terms  dimensionally, 


yfct)  = 
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m  lo2- 


SlH  \oj 
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(A. 3. 8) 


which  corresponds  (as  it  should)  to  the  solution  obtained  by 
Hussain  and  Reynolds  (29)  where  the  bending  and  damping  force 
terms  were  not  even  considered.   This  simplification  of  the 
solution  verifies  that  the  central  portion  of  the  ribbon 
essentially  behaves  as  a  massless  and  completely  flexible  vibra- 
ting string,  satisfying  the  equation, 


\rr\ 


J,=  K^   ^IBslocot 
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with  conditions, 


y(o,t)-   y(l.-t)-  0 


D.   Attainment  of  the  proper  mode  and  amplitude  of  oscillation 

To  this  point,  no  attempt  has  been  made  to  examine  the 
character  of  the  resulting  solution  for  ribbon  displacement  (i.e., 
points  of  resonance,  number  of  nodes,  effect  of  ribbon  tension, 
etc.)  or  to  answer  the  second  equestion  posed  in  section  4.2 
regarding  the  magnitude  of  the  force  required  to  drive  the  ribbon. 

Examination  of  equation  (A. 3.8)  reveals  that  resonance 
would  occur  when 


u,^  l  -_  Trn  n=    1,2,3,... 


K 
or         a;  -  jm  IK_ 

where  go   is  designated  the  "natural"  frequency  of  oscillation. 

Since  it  is  desired  that  the  ribbon  oscillate  in  a  single 
loop  standing  wave  as  visually  verified  in  use  in  other  experiments 
(7,  29),  parameters  affecting  the  number  of  maxima/minima  must  be 
adjusted  until  this  configuration  is  attained.   By  maximizing  the 
equation  (A.  3.8)  (and  realizing  that  ?  /I  6.   1  ),  it  can  be  shown 
that  multiple  loops  form  only  if 
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WK^-  1  *2irn 


n  =  1,  2,  3,  . .. 


(The  n  =  0  case  is  just  the  desired  single  loop.)   Apparently, 
then,  multiple  loops  should  not  occur  until  after  having  passed 
through  the  second  resonant  condition.   After  that  point,  for 
2rr<  °°  ]-77  1  ^  4tt   three  loops  (or  four  nodal  points)  character- 


ize  the  motion,  for  4?r<  °°  l  Lf  1 <  »TT  five  loops  should 
theoretically  be  observed,  and  so  on.   Figure  4.3  shows  the  first 
two  modes  of  oscillation. 

With  the  formula  for  \A   (  2  ,  i  ) ,  the  force  required  to  drive 
the  ribbon  at  a  certain  frequency  for  a  desired  maximum  amplitude 
can  easily  be  determined,  keeping  in  mind  that  u)  j-?/"  1  <  2TT 
for  a  single-loop  oscillation.   This  amplitude,  occurring  at  the 
ribbon  midpoint  can  be  found  from  equation  (4.3.8)  to  be: 
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the  magnetic  flux  intensity,  £>  ,  can  be  determined. 


g  =        lj^kx  uo    rv\         _ _  _    9^53  gauss 


I   tan(«,lffljtan(^ 


Note  for  this  example,  W|  — ]  =  .9064  <£  Zv   and  so  only  a  single 
loop  should  be  expected. 

Just  as  the  maximum  and  minimum  allowable  temperatures  for 
the  existing  testing  configuration  were  computed  in  section 
3.2A,  limitations  are  imposed  on  the  range  of  frequencies  required 
According  to  Wazzan,  Okamura,  and  Smith  (reference  (2),  table  2 
and  figure  30) ,  the  nondimensional  maximum  frequencies  of 
unstable  disturbances  for  the  previously  specified  temperature 
interval  (0  ^  Tw  "  ~Q>  -  22.5°F,  7^  =  60°F)  are 


-3 
Unheated  plate        uj     =  .  23  5  V /0 

Tw-T„-  22.5°F    3    =  .235  x/o-  + 


where  co     =   ^  ^ 

woo 

Again  using  as  the  limiting  tunnel  speeds,  \j    =  40  fps  and 


oo 


IV  5  fps, 
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maximum  frequency  =  5980  Ha 
minimum  frequency  =9.4  \At 

4.4   Implementation  of  the  Method  for  Vibrating  the  Ribbon- 

Since  other  experiments  which  used  a  vibrating  ribbon  to 
perturb  the  flow  were  conducted  in  a  wind  tunnel,  very  little 
difficulty  was  encountered  in  putting  a  given  amount  of  tension 
on  the  ribbon  -  using  rubber  bands  and  weights.   In  this  design 
where  problems  with  sealing  against  water  leakage  from  the  tunnel 
were  encountered,  special  provisions  had  to  be  made  to  counter 
these  problems  while  still  permitting  adjustment  of  the  ribbon 
tension.   Rather  than  using  suspended  weights,  the  design  shown 
in  figure  4.4  was  used  -  a  glorified  spring-loaded  fish  scale! 
Note  that  while  tightening  of  the  bolt  in  the  plexiglas  frame 
to  compress  the  packing  (which  sealed  around  the  ribbon)  will 
seal  the  hole,  it  will  also  tend  to  give  erroneous  readings  on 
the  calibrated  tensioning  device  if  adjustment  is  attempted  while 
the  packing  is  compressed. 

This  is  not  felt  to  be  too  great  a  drawback,  since  tension- 
ing could  be  made  prior  to  sealing  tightly.   Once  the  tension  is 
set,  it  should  not  need  to  be  readjusted.   In  a  similar  line  of 
thought,  the  edges  of  the  plate  would  have  to  be  sufficiently 
rounded  at  the  points  where  the  ribbon  would  pass  so  as  not  to 
excessively  bend  the  latter. 
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To  eliminate  possible  electrical  current  paths  other  than 
through  the  wire  ribbon  (a  potential  shock  hazard) ,  the  3/4" 
outer  cylindrical  part  of  the  device  was  threaded  into  a  teflon 
insert  which  had  been  press-fit  into  the  aluminum  supporting  block. 
This  cylinder,  in  turn,  was  encased  in  5/8"  x  3/32"  clear  tygon 
tubing,  and  tensioning  was  accomplished  by  a  teflon  "nut"  at  the 
top  of  the  cylinder.   A  previously  noted  precaution,  specifically 
anodizing  the  plate  to  create  a  dielectic  surface,  eliminated  any 
path  from  the  wire  to  the  plate  (see  appendix  A-2).   Since  the 
packing  device  is  imbedded  in  the  nonconducting  plexiglas  frame, 
no  current  leakage  should  exist  to  the  water  tunnel. 

The  inability  to  easily  vary  the  strength  of  the  magnetic 
field  was  sacrificed  for  the  elimination  of  the  power  source 
and  temperature  control  devices  required  for  the  proper  operation 
of  electro-magnets  (see  (29)  for  example).   Hence,  two  sixty  pound 
pull,  ALNICO,  permanent  magnets  were  chosen.   Since  it  was  realized 
that  the  magnetic  flux  density  would  drop  rapidly  away  from  the 
magnet  face,  a  recess  was  milled  in  the  3/8"  thick  flat  plate  to 
place  the  magnets  to  within  1/16"  of  the  working  side  of  the  plate. 

With  other  design  considerations  in  mind,  it  was  determined 
that  the  closest  that  the  ribbon  could  be  placed  to  the  leading 
edge  was  2.5".   These  ribbons  were  made  from  spring  phosphor 
bronze  due  to  that  material's  high  fatigue  resistance,  tensile 
strength,  and  electrical  conductivity.   They  were  initially  .1" 
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wide  and  .0035"  thick,  and  stretched  across  the  width  of  the  plate, 
9".   It  is  anticipated  that  perhaps  the  width,  length  of  vibrating 
span,  and  distance  from  the  plate  might  need  to  be  adjusted  to 
produce  a  measureable  disturbance,  although  Schubauer  and 
Skramstad  (7)  indicated  that  for  operation  in  a  wind  tunnel,  none 
of  these  factors  significantly  affected  the  performance.   Naturally, 
the  further  into  the  boundary  layer  the  ribbon  remained,  the  less 
flow-  and  self-induced  flutter  should  be  expected.   Although  the 
amplitude  of  these  disturbances  could  be  analytically  calculated 
from  the  previously  discussed  model,  it  is  thought  that  the  ampli- 
tude constituting  an  "infinitesimal"  disturbance  could  only  be 
determined  by  experimental  measurement. 
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APPENDIX  A 
ANODIC  COATING  ON  THE  ALUMINUM  BASE  PLATE 

A.  1   Formation  of  Anodic  Coating 

The  anodizing  process  used  was  the  ALCOA  Alumilite  coating  226, 
standardized  as  2  mil  thickness.   The  properties  of  the  dense,  ad- 
herent anodic  oxide  coating  to  be  discussed  in  section  A. 2  are  a 
function  of  the  electrolyte  solution  composition  (a  12%  sulfuric  - 
1%  oxalic  acid  solution),  bath  temperature,  composition  of  the 
alloy,  the  current  density,  and  the  elapsed  time.   As  pointed  out 
by  Wernick  and  Pinner  (20) ,  the  actual  mechanism  of  anodic  oxida- 
tion is  very  complex  and  controversial,  and  dependent  on  the  types 
of  electrolytes  used  to  determine  the  composition  of  the  coat- 
ings.  In  simple  terms,  in  this  electrochemical  process  an  electric 
current  is  passed  through  an  aqueous  electrolyte  solution  in  which 
the  aluminum  piece  to  be  anodized  acts  as  the  anode  and  the  lead 
walls  of  the  container  act  as  the  cathode.   The  negatively  charged 
anions  consist  in  part  of  oxygen  ions  which  react  chemically 
with  the  aluminum  to  form  the  anhydrous  aluminum  oxide  in  accordance 
with  the  formula: 


2  Al  +  ^  0  —*A\Ot     ♦  Energy 


Since  the  reaction  products  are  sparingly  soluble  in  the  electro- 
lyte, the  film  growth  on  the  surface  of  the  anode  is  accompanied 
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by  its  dissolution  as  well.   It  should  be  noted  that  the  oxalic 
acid  is  added  to  the  electrolytic  solution  to  decrease  the  solu- 
bility of  the  coating  serving  to  increase  its  hardness  and 
limiting  thickness.   This  dissolution  forms  pores  in  the  coating 
which  are  wide  enough  to  allow  continued  access  of  the  current 
to  the  metal  (specifically,  the  pore  diameter  is  on  the  order  of 
330°  A  or  3.937  x  10~9  inches  (33)).   The  inward  film  growth, 
retarded  by  the  increased  electrical  resistance  caused  by  the 
thickening  film,  continues  until  the  rate  of  growth  and  dissolu- 
tion are  equal,  at  which  point  the  film  thickness  remains  constant 
The  coating  will  also  contain  small  percentages  of  sulfate  and 
oxalate.   In  addition  to  the  primary  reaction,  a  number  of 
secondary  reactions  also  occur  at  the  anode  such  as  the  reduction 
of  oxalic  acid: 

(CooHl  +  0 2  C0L  +  HO 

and  the  electrolysis  of  water  (40) .  Gas  analysis  at  the  anode 
has  indicated  the  presence  of  carbon  dioxide,  carbon  monoxide, 
oxygen,  hydrogen,  and  nitrogen  with  proportions  depending  upon 
the  current  density  (20). 

After  plugging  desired  holes  with  rubber  or  plastic  plugs 
to  restrict  the  area  to  be  anodized,  the  oiece  was  then  put  in 
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the  electrolyte,  ensuring  positive  electrical  contact,  and  the 

current  turned  on  raising  slowly  2-5  volts  per  minute  to  the 

2 
prescribed  current  density  (in  this  case  36  amp/ft  ).   Due  to 

the  discharge  of  energy  in  this  reaction  as  indicated  above,  the 

electrolyte  had  to  be  agitated  and  refrigerated  (8-10°C)  since 

insufficient  agitation  would  cause  locally  higher  electrolyte 

temperatures,  which,  in  turn,  could  cause  burning.   The  thinner 

the  sheet,  the  more  agitation  was  required  to  prevent  burning. 

A. 2  Effect  of  Coating  on  Plate  Performance 

An  assessment  of  tho  various  properties  of  the  anodized 
coating  and  how  the>  would  affect  the  performance  of  the  plate 
was  necessary.   A  summary  of  that  study  is  presented  below. 
Hardness  and  abrasion  resistance 

An  assessment  of  the  wearing  quality  as  measured  by  con- 
ventional indentation  tests  are  not  indicative  of  the  material's 
hardness,  since  such  test  points  load  the  coating,  forcing  it 
into  the  much  softer  Al  substrate  producing  too  deep  an  impression 
for  reliable  results.   For  example,  repeated  sharp  blows  from 
a  hammer  can  dent  and  slightly  chip  the  coating,  owing  to  the 
thinness  of  the  coating  and  the  softness  of  the  basis  metal. 
However,  we  have  already  specified  that  the  plate  would  not  be 
subjected  to  localized  loading,  so  we  turn  to  the  results  of 
two  abrasion  tests  -  resistance  to  wear  by  rubbing  and  resistance 
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to  corrosion  from  a  grit  blast.   In  such  tests  the  resistance  to 
abrasion  of  the  hard  anodize  coat  was  found  to  be  about  twice 
that  of  the  regular  coating  (169  gm/mil  of  coating  to  89  gm/mil 
for  std  coating,  13  times  that  of  untreated  Al ,  and  comparable 
to  case-hardened  steel  and  chromium  plate.   For  rubbing  wear, 
6061  T6  hard  coated  Al  performed  about  like  case  hardened  steel 
and  much  better  than  carbon  steel  or  type  304  stainless  steel. 
Resistance  to  corrosion 

Even  though  the  pores  formed  in  the  hard  anodizing  process 
are  not  sealed,  since  that  would  reduce  the  abrasion  resistance, 
the  coating  provides  excellent  resistance  to  atmospheric  and 
saltwater  corrosion.   Compared  with  conventional  coatings,  which 
are  usually  expected  to  withstand  a  few  hundred  hours  of  salt 
spray  exposure  without  failure,  6061  alloy  with  a  2  mil  hard 
coat  withstood  10  months  exposure  as  per  ASTM  B117-49  with  no 
sign  of  failure  (15).   Experiments  on  the  effect  of  rain  erosion 
on  the  leading  edges  of  high  speed  airfoils,  conducted  by  Gillig 
show  that  failures  depend  primarily  on  the  alloy  which  is  hard 
coated.   The  best  results  were  obtained  with  24S  Alclad,  61S, 
and  24S,  in  that  order  (14). 
Fatigue  and  tensile  strength 

The  anodizing  process  causes  a  slight  decrease  in  the  ten- 
sile strength  (for  61S  T6  with  a  2  mil  coating  a  reduction  from 
47,700  to  47,000  psi)  believed  due  to  the  stress  concentrations 
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at  the  microcracks  in  the  coating.   In  general,  since  anodic 
coatings  in  excess  of  .5  mil  tend  to  reduce  the  fatigue  strength 
of  wrought  Al,  hard  anodic  coatings  are  generally  not  used  when 
the  fatigue  strength  is  a  factor;  but  for  this  application,  it  is 
not  a  factor.   The  slight  loss  in  ductility  is  not  considered 
important. 
Adherence 

Due  to  the  nature  of  the  process ,  the  coating  becomes  an 
integral  part  of  the  piece.   The  only  question  of  adherence  would 
be  at  the  leading  edge  as  discussed  above.   Although  severe  bend- 
ing of  the  piece  should  not  be  experienced,  tests  show  that  a 
right  angle  bend  over  a  3/4  inch  diameter  shows  that  a  2  mil 
coating  spalls  off  on  the  compression  side  of  the  bend  and  checks 
finely  on  the  tension  side,  but  the  piece  can  be  bent  180  degrees 
on  a  liberal  radius  and  straightened  out  again  without  spalling 
off  of  coating.   The  adhesion  of  the  oxide  coating  is  normally 
much  better  than  that  of  electrodeposits. 
Electrical  properties 

Since  the  wire  ribbon  of  the  perturbing  mechanism  could 
come  into  contact  with  the  plate,  it  would  be  advantageous  to  have 
the  latter  insulated.   For  the  nonconducting,  unsealed  anodic 
coating  (2  mil),  the  breakdown  voltage  is  950  volts.   Although 
sealing  in  boiling  water  of  the  coating  increases  the  breakdown 
voltage  for  2  mil  coating  to  1,200  volts,  it  is  not  considered 
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necessary  to  obtain  the  additional  resistance  at  the  expense  of 
the  abrasive  resistance  (although  the  plate  would  still  be  well 
within  the  limits  required  for  this  application  in  either  the 
sealed  or  unsealed  configuration). 
Effect  on  surface  finish 

Due  to  the  excellent  throwing  power  of  the  process,  uniform 
coatings  may  be  applied  to  irregular-shaped  objects  but  varying 
slightly  for  a  closed  end  hole  (i.e.,  the  thickness  of  the  coating 
in  the  unplugged  blind  bolt  holes  in  the  plate  could  not  be 
expected  to  be  uniform) .   In  general  the  surface  roughness  is  slightly 
increased  by  hard  anodizing.   For  example,  for  7075  T6  with  a 
coating  thickness  of  1.2  mil,  an  initial  surface  smoothness  of  5-7 
microinches  (RMS)  will  have  a  smoothness  of  10-20  microinches 
after  the  treatment.   If  a  finer  finish  were  desired,  the  plate 
could  have  been  made  oversized  and  then  lapped  or  honed  down  to 
a  smooth  surface  on  the  order  of  2-3  microinches.   Since  the  plate 
was  machined,  then  lapped  prior  to  treatment,  it  is  felt  that 
the  surface  smoothness  is  well  within  the  limits  at  which  it 
will  contribute  to  turbulence  generation. 
Heat  resistance  and  thermal  expansion 

In  general,  due  to  the  differences  of  coefficients  of  thermal 
expansion  between  the  anodic  coating  and  the  aluminum  substrate, 
checking  of  the  film  occurred  when  the  temperature  exceeded  200°F. 
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The  coating  itself  acts  as  an  insulator  -  a  disadvantage  when 
trying  to  conduct  heat  to  the  surface  of  the  plate  to  be  con- 

vected  into  the  boundary  layer  -  having  a  thermal  conductivity 

2 
of  7.17  Btu/(hr)  ft  F/in.  (41).   The  coating  does,  however, 

have  a  good  resistance  to  flame  impingement  or  short  term  high 

temperature  exposures.   When  heat  transfer  is  by  radiation 

emissivities  of  90-95%  are  to  be  expected.   Due  to  the  highly 

refractive  nature  of  the  coating,  the  interior  will  melt  before 

the  coating  collapses. 


APPENDIX  B 

ANALYTICAL  SOLUTION  OF  HEAT  CONDUCTION  IN  A  COMPOSITE 
MATERIAL  WITH  HEAT  GENERATION  IN  ONE  REGION 

Heat  is  generated  by  a  finite  series  of  rectangular  heat 
sources  in  one  region  of  a  two-region,  composite  rectangular  slab 
which  is  insulated  on  three  sides  but  which  transfers  heat  by 
convection  to  a  laminar  boundary  layer  on  the  remaining  side. 
It  is  required  that  the  temperature  distribution  on  the  plate 
surface  adjacent  to  the  flow  be  found  so  that  the  influence  of 
spacing  and  strength  of  the  sources  can  be  evaluated. 
B . 1  Determination  of  Source  (Neumann)  Function 

Following  the  notation  of  the  mathematical  model  shown  in 
figure  B.l,  the  equation  to  be  solved  with  accompanying  boundary 
conditions  can  be  formulated: 

(where  the  point  source  discontinuities  are  in  region  a  ),  with 
boundary  conditions 

K*^(V)*    ^(x,o)  M^o)-  Mc(X;o) 
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Note  in  the  specification  of  boundary  conditions  that  the  Neumann 
function (s)  is  (are)  derived  with  homogeneous  boundary  conditions 
(see  (A2)  for  elucidation). 

Bulavin  and  Kashcheev  (A3)  have  solved  the  transient  heat 
conduction  problem  for  multilayer  symmetrical  bodies  with  heat 
generation  within  the  solid  by  the  method  of  separation  of  variables, 
and  by  construction  of  orthogonal  expansion  of  functions  over  multi- 
layer regions  for  the  case  of  perfect  thermal  contact  between 
the  surfaces.   Similarly,  for  the  steady  problem  here,  it  is 
assumed  that  (\|»(  X  ,  M  )  can  be  expanded  in  an  orthogonal 
eigenf unction  expansion, 

OO  CO 


NjCx^-x^OO     )     ^nNm(x'^')NJn(3)cos 


"2C 

<  -  JL 

on      2  (B.1.3) 


nnn 


=   1  YY\f   O 


where   f     are  eigenf unctions  satisfying  the  following  differential 

•J 

equation 


Now  considering  the  orthogonality  property  of  the  eigenf unctions 
Y,    over  the  entire  rectangle,  i.e.,  over  the   I   layers,  where 
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rtn 


K; 
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Tl 


J 


w 


lJnCy)  lj-'^)^  =  (constant)  <Snn<    (B>1*5) 


with    !(".   being  the  steady,  discontinuous  weighting  function 
(thermal  conductivity)  for  layer   i  which  is  unchanged  for  either 
the  assumption  of  perfect  thermal  contact  or  linear  contact 
resistance  at  the  interface  or  for  boundary  conditions  of  the  first, 
second,  or  third  kinds  (44). 


ra. 


J 


r^ 


N-.  Cos-       ■  Kj  \-n  du  ax   = 


* 


rWi+i 


Kj  |\Jj   Y:n  Cos  L^L  ciuc/x 


5j 


ij,  «, 


(B.1.6) 


where         M  -    )     < 


'J 


J 


--1 


^  d7 


(B.1.7) 


NmnC^*j#)  = 


2. 


••      jTl^ 


(B.1.8) 


k  Nj^oos^djdx 
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Applying  Integral  transform  (B.1.8)  to  equation  (B.l.l), 
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O        0 


K    J  a  'an  1 


■>3L  r<>   /    2- 


o    -c 


(B.1.9) 


but, 


^x2  a.  I  a>x  * 


or  after  applying  boundary  conditions  from  equation  (B.1.2) 
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ttii 


J  ^X2  #• 


TTrr\ 
3L 


)  J  Na  cos 


^^  dx 


(B.1.10) 


Similarly 


CcS 


ITvnx 


dx  =  -  (™)  J  N(c  c*  2l™L  d  x 


(B.l.ll) 


also 


^j1    'an      J      y    Rn    c)W 


1 


b       ,b 


+ 


Ill 


or  applying  equations  (B.1.2)  and  (B.1.4) 
du    'an   J        an    £u   J 


■+- 


Similarly 

-c  Y 


^(M*®-^  YMN.«A( 


Ycndr  Ycn(o)^(x,o) 

-Nc(X)0)d^a(o)-^j_°    YcnN0c/j 


(B.1.12) 


(B.1.13) 


Substituting  equations  (B.1.10)  through  (B.1.13)  into  equation 
(B.1.9), 

a  r^ 


aM 
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+  ^ 


^NX^s^cl^, 


o  -o 


f(^x-^^d,]+^{r 


-OU°)^ft$ 


J  J 


dYc 


:  dxf=  -_2_ 


-5T 


TTyv\X 


or  upon  simplification  with  equations  (B.1.2)  and  (B.1.8) 


•Y„op 
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NU^')  = 


xji')** 
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TT^ 


+■ 


(B.1.14) 


4m 


With  equations  (B.1.3)  and  (B.1.14)  thus  established,  the 
problem  is  essentially  solved;  all  that  remains  is  to  determine 
the  expressions  for    I'm  C  H  )  anc*   (2>«  ♦   Recalling  equation 


(B.1.4)  where 


t" 


v%(y)*£  V:Ps0 


it  can  be  shown  that  the  solution  takes  the  form 


V:P"Aj«<M)  +  6r%V 


where    (U-  (  t|  )  and    y'n^  ^  )  are  two  linear>ly  independent 
solutions,  which  in  rectangular  coordinates  take  the  form 


^j„G|)=    cos^jj 


Therefore, 


Yn(^j)  -  Ajn  5'm^ij  4-   BJ0  cos  ^u  (b.i.15) 
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Since    j .   ( U  )  must  satisfy  boundary  conditions  in  equation 
(B.1.2),  equation  (B.1.15)  can  be  substituted  into  these  boundary 
conditions  to  yield 


or 


and 


<^toS^b    -     Bar,5'°^b    =0 

6a„  -  6,    =  0 

^cnCo5^C  +   B-  S,°  ^C    =  ° 


^C  +    DCn  sin  6r 


&cn=   Ban=  A^cct^b 


(B.1.16) 


Acn  "    -rr-  Aan  (B.1.17) 

tsc 


CoU    b     +  y1   Cot  a   C    =    0  (B.1.18) 

Hence,   Yi'n  (  U  )  is  determined  in  terms  of  one  constant 

A 
(here  chosen  as   f»an  )  which  will  cancel  as  the  solution  progresses, 

and  6n  are  the  eigenvalues  determined  by  equation  (B.1.18). 
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From  equations    (B.1.15)    and   B.1.17) 


V      -     f\         Cos  ft^(tj-b) 


(B.1.19) 


en 


Now,  to  determine     from  equations  (B.1.7)  and  (B.1.19),  with 
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(B.1.20) 
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Combining  equations  (B.1.3),  (B.1.14),  (B.1.18),  (B.1.19) 
and  (B.1.20)  the  source  functions  are  completely  determined,  and 
may  be  written  as 


00     00  / 


U    -    A    \    \       J„„  CoSftrt(if-b)coSftnfu-b)5cnftnC  Cos  ^ Cob 


and  (B.1.21) 


nn»0  0=1 


nvrwi     ciTnf\y 


(^)V^][bkn^ct|asm^bj 


where  X        =1  m  4  0 

rr\r\  ' 


B.2  Determination  of  Temperature  Distribution 

With  the  source  functions  now  determined,  attention  is  directed 

toward  finding  the  temperature  distribution  on  the  surface  of  the 

rectangular  slab  resulting  from  a  series  of  rectangular  heating 

elements  of  strength  <S-   ,  thickness      ,  and  width  A   ,  lo- 
co. ^ 

cated  at  positions  (  ><•  ,  M'L  )  in  region  <3.  .   To  represent  these 
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sources,  a  combination  of  step  functions  is  used  so  that  the 
internal  source  distribution  is 


L 

H[r(^-4)]-H[r(^-4)]] 


(B.2.1) 


where  it  is  assumed  that  o.„   is  uniform  within  each  element, 
and  that  the  dimensions  on  all  elements  are  identical.   Using 
the  procedure  outlined  by  Duff  and  Naylor  (42),  with  the  kernel 
or  Neumann  source  function  now  defined  by  equation  (B.1.21), 
the  temperature  can  be  determined  to  within  an  additive  con- 
stant with  any  desired  internal  source  distribution  and  speci- 
fied surface  heat  flux  as  follows: 


9  =  T-T; 


ave. 


N;s.dA  + 
J  J  J 


M.ftJ! 


(B.2.2) 
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where       I  =  local  temperature 


=  average  temperature  over  the  region  =  —  \~T~dlh\ 


o«  =  internal  source  in  region  I 


t 


■     j    on  the  external  boundaries  of  region  | 

J    ~dr\  J 


subject  to  the  restriction  that  in  a  steady  state  configuration 
the  heat  generated  internally  must  be  transmitted  through  the 
boundaries. 


(B.2.3) 


If  this  condition  is  not  satisfied,  then  no  solution  can  exist. 

Recalling  the  discussion  of  the  model  to  be  used  in  this 
analysis  as  presented  in  the  introduction  to  this  appendix,  the 
following  boundary  conditions  are  stipulated: 

=    -    H_ 

(B.2.4) 


Pe(VC) 


T^ 


P^)  = 


x  0 


on  all  other  boundaries 


k,\4 


where  for  want  of  a  better  choice  of  a  constant  |-|  =  .332  K  Vr  [  °°)  AT 


v 


in  accordance  with  the  similarity,  Blasius-type  solution  for  a 
thermal  boundary  layer. 
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Incorporating  eqns.  (B.1.21),  (B.2.1),  and  (B.2.4)  into 
equation  (B.2.2),  the  temperature  distribution  in  the  plate 
would  be 


a_    4  aVX  Vco5pn(w-b)SLnafioCa)Spn(u.-b)  Sin  (5 
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4- 


A 


TTnAy 
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L(^fv, 


2.9  cos  TI^  sin  TL^  cos^  ( tj:  -  b)cc^n  (ij-k)sLn  ^| 


(B.2.5) 


a  prodigious  looking  equation  to  say  the  least! 

Furthermore,  to  insure  the  validity  of  a  solution,  the  con- 
dition specified  by  equation  (B.2.3)  must  be  satisfied.   Therefore, 
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JsjvpdA  =-jpc(x>-^)d( 


V^AE  =  ZUfzC 


or  A     AR--     7U/7  (B.2.6) 


After  attempting  to  determine  the  eigenvalues  with  the  charac- 
teristic constants,  it  becomes  apparent  that  the  eigenvalues  are 
approximately  (3n~ — r~    n  -   0,  1,  2,  3,  ...  .   Specifically, 
for  b-  .25",   C  =  .002",   K^=  117  Btu/hr  ft  F,  and    fy  -  .597 
Btu/hr  ft  F, 


61  »  3.999830   ~  ^iXll 

(3  =  7.999673   ~   -Ullfl 
V2  .25 


It  can  be  seen  that  the  real  eigenvalues  will  be  less  than  the 
approximate  ones  (a  difference  that  will  increase  as  n — *-co   ), 
but  owing  to  the  relative  rapidity  with  which  the  solutions  are 
expected  to  converge,  this  difference  is  not  expected  to  signi- 
ficantly affect  the  solution. 

Therefore,  specifying  that   Qn  ■  — r—  ,  the  source  functions 
of  equation  (B.1.21)  are  drastically  simplified,  becoming, 
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(B.2.7) 


where 


x--  i 


rr\n 


I. 


*oo  ~-   O 
Physically,  this  implies  that  the  lateral  heat  transfer  in  region  C 

is  insignificant  when  compared  with  that  in  region  3.  .   Consequently, 

the  slab  should  be  treated  as  a  single  component  material  and  the 

source  function,  or  kernel  N -»  be  used  with  all  sources  and 

a. 

boundary  conditions  in  the  integral  equation  (B.2.2). 
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This  equation  can  be  simplified  somewhat  by  evaluating  some  of 
the  series.   From  Jolley  (45) 
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Making  the  appropriate  substitutions,  the  temperature  distribu- 
tion is  finally  obtained. 
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B.  3  Numerical  Examples  of  Temperature  Pis trlbu^t ions  for  a  One 
Material  Mod el  for  Specified  Heat  Source  Distrlbut ions 

As  an  example  of  the  solution  to  equation  (3.2.9),  the  required 
heat  flux  was  assumed  to  be  provided  by  twelve  sources.   An  initial 
indication  of  the  strengths  and  positions  of  these  sources  is  ob- 
tained from  equation  (3.2.14),  derived  in  section  3.2C,  for  the 
following  configurations: 

1.  Constant  strength  sources  with  variable  spacing 

2.  Variable  strength  sources  with  constant  spacing 
Pertinent  parametric  values  applicable  to  both  are  found  in 
table  3.1.   Source  strengths  and  positions  in  the  plate  for  con- 
figurations (1)  and  (2)  above  are  found  in  tables  3.2  and  3.3, 
respectively. 

In  the  evaluation  of  equation  (B.2.8)  for  a  specific  example, 
the  Fresnel  integral  is  encountered,  i.e.,  an  integral  of  the 
form 
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which  must  be   approximated  as   follows    (see  Jolley    (45)): 
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In  practice  the  first  formulation  is  used  for  rr\  £  3  and  the 
second  for  all  higher  values  of  fY]    . 

To  obtain  the  desired  accuracy  in  the  temperature  distrubu- 
tion,   K  9  of  equation  (B.2.8)  is  summed  on  the  computer  so 
that  the  individual  infinite  series  are  converged  to  an  accuracy 
within  1.01.   The  FORTRAN  V  program  used  is  shown  in  figure  B.2. 

Surface  plate  temperatures  resulting  from  the  example  con- 
figurations (1)  and  (2)  are  shown  in  figures  3.2  and  3.3,  respectively, 
with  characteristic  details  of  the  latter  appearing  in  figure  3.4. 
A  discussion  of  the  ramifications  of  these  examples  and  figures 
is  presented  in  section  A. 2D. 
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APPENDIX  C 
ASYMPTOTIC  SOLUTION  TO  VIBRATING  RIBBON 

C . 1  Formulation  of  the  Problem 

A  wire  ribbon  of  length  1,  width  C~  ,  and  thickness  K  ,  is 
stretched  parallel  to  a  flat  plate  at  a  close  distance  and 
constrained  at  each  end  with  built-in  end  conditions  (i.e.,  the 
slope  and  deflection  at  the  end  points  are  zero).   It  is  forced 
to  vibrate  electromagnetically  in  a  viscous  medium  by  passing  an 
alternating  electric  current  through  the  ribbon  of  frequency  (jl) 
in  the  presence  of  a  constant  magnetic  field.   To  determine  the 
effect  of  small  viscous  damping  and  bending  on  the  motion  of  the 
ribbon,  the  above  model  is  formulated  and  solved  with  the  desired 
effects  being  just  perturbations  on  the  main  forced  oscillation. 

Consistent  with  the  notation  in  the  mathematical  model 
shown  in  figure  C.l,  the  defining  equation  is 


(C.l.l) 

U y~ J  ^ y *  »> y '       I- J- J  <■ 1 

bending    inertial       viscous     tension     forcing 
force      force         drag        force       function 


Nondimensionalizing  with 

r=u>t       f=X       *SH 
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equation  C.l.l  becomes 


+  ^ 


i£ 

br7 


where 


af 


=  "KSinT 


(C.1.2) 


X-   IB1 


60' 


*2-  fr-^r 


60. 


oo. 


1       nn 


i-  t'cftm*' 


with  boundary  conditions; 

spatially:    -|£- ^  f1  =  0   at  ^  =  °,  1 

temporally:  the  requirement  that  the  viscous  damping  be  a 
maximum  when  the  ribbon  velocity  is  a  maximum 
and  zero  when  the  ribbon  is  locally  at  rest. 

Note  that  the  resulting  equation  is  a  fourth  order,  nonlinear, 
partial  differential  equation.   Examination  of  this  equation 
suggests  two  possible  ways  of  reaching  an  approximate  solution: 
1)  a  regular  perturbation  in  A  ,  solving  the  fourth  order 
equation  exactly,  or  2)  a  boundary  layer  type  solution,  with 
interior  and  exterior  expansions  matched  in  terms  of  an  inter- 
mediate variable.   The  interior  solution  would  be  expanded  in 
a  regular  perturbation  expansion  in  6  . 
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The  first  approach  leads  to  such  a  complex  zeroth  order  solution 
(no  damping)  with  the  end  effects  incorporated  within  it,  that 
the  higher  order  corrections  (where  the  effect  of  the  damping  is 
considered)  prove  to  be  too  cumbersome  for  our  purposes  here. 
Specifically,  since  the  behavior  of  the  wire  ribbon  near  its 
center  -  away  from  the  end  points  -  is  the  quantity  of  interest, 
the  second  method  suggested  above  is  used. 

C . 2  Asymptotic  Solution  for  Interior  Region 

With  the  viscous  damping,  it  is  expected  that  perhaps  there 
might  be  more  than  one  time  scale  associated  with  the  problem,* 


that  is,  one  associated  with  the  damping  (  T  =  £>T  )  and  one 


with  the  oscillation  of  the  wire  (  T  =T  (  1  ■+-  6  t02j  )  where  k02 
is  a  frequency  shift  caused  by  the  damping. 

Following  Cole  (46),  it  is  assumed  that  for  the  interior 
portion  of  the  ribbon,  the  dependent  variable  of  the  defining 
equation  can  be  expanded  as 

'        V  v         (C.2.1) 

where  \)  (t)  \)  (£)...  are  some  functions  of  £  to  be  determined 
later.   By  substituting  equation  (C.2.1)  into  equation  (C.1.2), 
the  ordered  set  of  equations  is  powers  of  €.  Is  obtained. 


132 


^-     +  J&H    +^h 


r 


dr: 


ar 


ak 


7? 


=  tSinT  (C.2.2) 


(Note  that  the  time  used  in  the  forcing  function  is  that  asso- 


ciated with  the  oscillation  of  the  wire.) 

r 


61.   -^1 


ar 


_  -< 


_ ak  if  i .  1 


o    if 


v>, 


(C.2.3) 


Further,  for  small  damping,  define 


(C.2.4) 


so  by  substituting  equation  (C.2.4)  into  equation  (C.2.2)  (for 
the  zeroth  order  solution  in  £  )  an  ordered  set  of  equations 
in  b    is  obtained. 


* 


i  x2- 


6°:    -  fL&i  +  *  Qfi*      =  ^  sua  r* 
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By  considering  the  zeroth  through  the  second  order  equations, 
utilizing  the  fact  (determined  later  by  matching  the  interior 
and  exterior  solutions  -  see  equations  (C.3.7)  and  (C.3.15), 
matching  conditions  (1)  and  (5)1  that  the  zeroth  order  solution 
in  £  must  satisfy  homogeneous  boundary  conditions  at  E  a  0,  1, 
and  requiring  that  the  solution  be  uniformly  valid  in  T      ,  it 
can  be  found  that  Po  through  C  ^  are  not  functions  of  T 
and  that   60.=  0.   This  means  that   T  =  T  ;  in  other  words, 
there  is  only  one  time  associated  with  the  problem.   In 
retrospect,  this  seems  reasonable  for; since  the  wire  ribbon 
is  driven  externally,  its  frequency  should  be  unaffected  by 
the  damping. 

Therefore,  equations  (C.2.1)  and  (C.2.4)  can  be  redefined 
and  combined  to  obtain  (to  first  order  in  (2>  ) 


+  .  .  . 


(C.2.6) 


Equations  (C.2.5)  for  zeroth  order  in  £  can  now  be  rewritten  as 


13A 


?; 


t: 


-  1&    +  <* 

>r2 


ar2 


o< 


ar 


ad*& 


dr: 


=    "*  sCn  r 


2)T 


ar 


(C.2.7) 


It  follows,  then,  that  the  zeroth  order  solution  in  &     and  6  is 


? 


i»(f,T)=    Sinr[A0Sm^|  4-    60Co5c<|    -X^J 


4  Cos  r[C0  sin  c^    f  D0  Coso( 


f] 


(C.2.8) 


C . 3  Asymptotic  Solution  for  the  Exterior  Regions 

and  Matching  with  That  for  the  Interior  Region 

To  determine  the  exterior  expansions,  it  is  assumed  that  the 
bending  and  the  tension  terms  are  of  the  same  order  of  magni- 
tude.  Stretching  the  E  variable  in  the  boundary  layer  area 
near   £  »  0  so  that 

the  governing  equation  (C.1.2)  becomes, 


14-4 


"  M      -  c1-  d'( 


~  A 

3f 


jLL^  +  B  df 


»?* 


^r 


^r 


rf 


7)T 


=  ^  Sin  t 
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By  equating  the  coefficients  of  the  hending  and  tension  terms, 
the  exponent,  fl  ,  is  found  to  be    n  =  -  1/2.   Therefore,  the 


boundary  layer  coordinate  near  ~f  =   0  is 


(C.3.1) 


Further  expanding 


I  (f;r5£)=^0(€)fo(f)r)  y.C^^p)^)^!^ 


■*-.. 


the  governing  equation  near   c  =  0  b 


ecomes 


4/< 


w-(a^;/ 


£IU...\ -i^jlLx  „  si, 


'^ 


f  .  • 


brz 


dT 


ty3 

dr 


a -ft 


dr 


^ 


^r 


+  .  •  •  \  ~    Ysinr 


To  lowest  order,  two  choices  are  available  for  the  form  of 


or 


l.  /-£- — > 


2. 


/' 


oo 


as   £  — ►oo 


(C.3.2) 


1 


Assuming  (1)  is  valid  (until  contradicted  or  confirmed) 


*%      A  . 


=  O 


(C.3.3) 


a? 


7 


By  satisfying  the  boundary  conditions  at   F  =  0,  and  assuming 
that  the  solution  of  the  form  £       will  not  match  with  the 
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interior  solution  as 


-^-  CO 


I-  A0cr)l^-Ue   ] 


(C.3.4) 


Proceeding  formally,  an  intermediate  variable  is  chosen 
suitable  for  matching  near  £  ■  0;  i.e.,  it  is  desired  that 


00  while 


? 


as 


Defining  ~fc       =  — C— v    where   T 


is  fixed  and  for 


e— >o 


1 


O 


as  £-*0    (C.3.5) 


In  preparation  for  matching  the  interior  and  exterior 
solution,  equations  (C.2.6)  and  (C.3.4),  the  former  is  expanded 
in  a  Taylor  series  expansion  about  Y    =0. 

+  . ..]  +  X-^  [&".  lo,r)+^%\  (o,t)-k.  .]  +  ... 
+  ^l(«)[K0(olr)+*K,(P,T)+...]  t  v^ylX^r) 


+-  pjc  co,t^+...J  +■  *  •  • 
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where  primes  indicate  differentiation  with  respect  to  spatial 
coordinate,   v 

Matching  the  interior  and  exterior  expansions  in  terms  of 
the  intermediate  variable  IE   near   E  =  0 

h        T 


1 


...]  +-...-hv)i(fe)[KD(o)r)+  sK^o/r)  ■»-...] 


r 
-1+e     J  yM,(6)  ^(^  T)-  ..  ^  =  0    (C.3.6) 


Matching  is  possible  only  if: 


1.    ^(O.^-r^^..-  0 

2.   A(€)  =  fi" 


(C.3.7) 


4-       A0Cr)  =  -[<o^r)-^K.(0^i-...] 

Note  that  since  M0(Am   lHP  ,  the  assumption  that  >^-  — >  co  is 
justified.   Note  further  that  the  term  of  0(1)  in  \       cannot  b< 


138 


matched  except  by  a  term  from  r\       ,  implying  that 

V  O)  -  |iT  (c.3.8) 


Since  . ' v  ;^>.co  £•— *o>  tne  zeroth  order  solution  in  Q   in  the 
first  order  solution  in  £  can  be  found  (see  equation  (C.2.3)). 

Following  the  same  procedure  at  the  other  end  of  the  wire, 
the  stretched  coordinate  is  defined  as 


?"""  ~  ~^~~^ 


-\- 


and  the  expanded  deflection  -f-  becomes 


+ 

(C.3.10) 


Putting  this  in  the  governing  equation  (C.1.2),  and  applying 
boundary  conditions, 


C- At(^^M-cf 


(C.3.11) 


Now  defining  an  intermediate  variable  so  that  near  V.     ■  1 
as  <f  — •"  1   ,  then  T  — ►  "Mas   £  — »-  O 

*  -  (F-0 

v-^  (c,3-12) 
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where  jr   is  constant  and  as 


o~($)—^0 


crfe) 


Therefore, 


pl+  ^o-(6)  —I 


(C.3.13) 


+ 


—  OO 


as  £ ^0  and  C7"(6r)->0  with  "?r  fixed. 

As  before,  matching  of  the  interior  and  exterior  solutions 
in  terms  of  the  intermediate  variable  is  accomplished  near  £  =  1 


+  1-e 


2^>r  -, 


.  .A  =  0    (C.3.14) 


Matching  is  possible  only  if: 


5-  £0(1/0+^, 0/r>...-  o 

6-  A0Cr)=  ^'(^r)  ^^iC,^-.-    (C3.15) 

7.    Afcr)=  Ko(1,rH^Kl0lr)4.... 
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As  noted  previously  when  showing  that  only  one  time  scale 
is  associated  with  the  problem,  the  zeroth  order  solution  in  £ 
must  satisfy  the  homogeneous  end  conditions  (matching  conditions 
(1)  and  (5)). 

With  the  main  contribution  to  the  solution  for  the  vibrat- 
ing ribbon,   PQ  ,  now  specified  by  equation  (C.2.8),  first 
order  corrections  due  to  the  influences  of  viscous  damping 
(represented  by  Q\     )  and  bending  (represented  by  K0  )  can 
be  evaluated  by  the  solution  of  equation  (C.2.7),  and  equations 
(C.3.8),  (C.2.3),  and  (C.2.4),  respectively,  with  matching  and 
boundary  conditions  determining  the  constants. 

C . A  Utilization  of  Matching  and  Boundary 

Conditions  to  Determine  the  Complete  Solution 

Utilizing  the  information  derived  in  the  preceding  two 

sections,  the  first  order  corrections  to  an  undamped,  completely 

flexible  oscillating  ribbon  due  to  viscous  damping  and  bending 

can  be  found. 

D  -  the  basic  solution 


Using  matching  conditions  (1)  and  (5)  for  equation  (C.2.8) 
and  realizing  that  each  term  in  an  asymptotic  series  should  in- 
dividually satisfy  the  boundary  conditions, 
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kl?^     (A0^^+^CoSo(V^)S^r  (C2.8) 


so   that 


go(°»r)^    [Bo-^]SinT"    +  D0^*T    a°  <C'4-1) 


£c  (t/r)  =   p05»n^  +  B0cos^-IUinr+ 


65ino(4-  Doc°sx 


cosr-o 


To 

be 

valid 

for  all 

r 

A0  = 

*       (1 
**- 

-  CoSo<) 
sin  o( 

60  = 

ce 

=     0 

D„ 

--   0 

I 


tan 


&< 


(C.4.2) 


so 


&(?,T> 


-   Y 


L   tan-y  Sm^-t-CoSo^-l 


SinT 


(C.4.3) 


Pt  -  viscous  damping  correction  to  basic  solution 

With  the  form  of  C0    now  determined,  equation  (C.2.7)  can 
be  solved  to  yield   C,  (  £  ,  T  ).   To  simplify  the  solution 
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the  nonlinear  damping    force,   being  an  even   function  in   T, 
is   expanded   in  a  Fourier  cosine  series. 


& 


br 


5r  ■ 


^r 


cos  nr 


where 


r\-\ 

~  -$-r  ta*n|-Sino<E  +  CoSo<£  -  \\  Cosr 


(C.4.4) 

COS  T 


b     =      Z   X 


-,2-r    rtt 


ta/\  y5in^+  Cos  o(  F  -  | 

ft* 


cosV  cos  rvrdf 


cos  r  cosnr  dr 


IT 


or 


b  =  ^L 


+ 


-.£,- 


tan^  Sin  c*jf  +.   CoSo^  -1 
Sm  C2-n)^_         Sin  (2+n) 


^^ 


n 


1T 


n=/2     (C.4.5) 


Therefore, 


2(2-n)  2(2tn) 


D  =  0  for  n  an  even  integer 


z  r  -.2. 

b,  a  |-  -t:  tan  2- sin  «S  +  Ccs*f  -1 


US  35 

b     -      -A- 


(C.4.6) 


Using  equations    (C.4.4)   and    (C.4.6),   equation    (C.2.7)   may  be 
written  as 
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-  &&L    +    ^L   t&L   -    -  k  (  Co^r  +    1   Cos  3' 


f 


3t 


5r  +  '    c«*7r-...) 


(C.4.7) 


Applying  the  matching  conditions  of  equations  (C.3.7)  and 
(C.3.15),  which  essentially  require  that  the  ribbon  satisfy 
homogeneous  conditions  at  its  end  points,  equation  (C.4.7)  yields 


+ 


5I45 


(tan^-<5in7^  +■  ^05  7o{f  -f)dos7r  -  ...J 


or 


where 


^3    '  875      \  tAo^-SiVwC    +-£o:>°\£-1/ 

K"      ft\  =      K1£\(t*^^.7*'l   +  to  7^-1 


(C.4.8) 
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\\    -   bending  correction  to  basic  solution 

With  t  (  £   >  T  )  now  evaluated  to  zeroth  order  in  £  , 
a  first  order  correction  due  to  bending  is  sought.   Since  the 
form  of  \).   (  G     )  is  determined  in  equation  (C.3.8),  equations 
(C.2.3)  and  (C.2.4)  become  (after  revision  necessitated  by 
the  presence  of  only  a  single  time  scale), 


V 


2  I         ^  IX 


a? 


2>T' 


aKc 


ar 


(C.4.9) 


From  zeroth  order  equation,       KQ    can  be   found   to   take   the 
following   form 


K0»    (A    Siho{[   +  B2CcSo<^)51nr+(C2Sin^f  +  D2CoSot^Cc5r     (C.4.10) 


From  matching   conditions    (3)    and    (4) 


K0(o,r)+-  aK^r)-*-,., 


so  that,  using  equations  (C.4.2)  and  (C.4.10),  for  zeroth  order 
in  6 

^  (ot)  =   ^  AaSinT=  -B^inT  -  U  CoS  T  -  -  K0(o,t) 
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Similarly,  using  matching  conditions  (6)  and  (7), 


|Xl)v)^_|(l)r)+...   =    Ko(lr)  +.  ^  K,(1,t) 


* 


+  . 


so  for  zeroth  order  in 


^•(1;T^  =  (A0*CoS<*-  QQ*Sif)*<)sinT  '     (A25Lino(  + 


B2cW)^nT  +  (Cz5Lnc^+ D^coso()co^r  =   Ko0>T) 


To  be  valid  for  all  time, 


A2  = 
B,= 

D2  = 


(C4.ll) 


0 

o 


Consequently,  with  equations  (C.4.10)  and  (C.4.11) 


^o(^T)---^t3.n(^)[tan(|)5Ln^^Cc5cv^5cnr 


(C.4.12) 
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